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An extensive investigation is made of the interior structure of self-similar accreting charged black 
holes. In this, the first of two papers, the black hole is assumed to accrete a charged, electrically 
conducting, relativistic baryonic fluid. The mass and charge of the black hole are generated self- 
consistently by the accreted material. The accreted baryonic fluid undergoes one of two possible 
fates: either it plunges directly to the spacelike singularity at zero radius, or else it drops through 
the Cauchy horizon. The baryons fall directly to the singularity if the conductivity either exceeds 
a certain continuum threshold Koo, or else equals one of an infinite spectrum k„ of discrete values. 
Between the discrete values k„ , the solution is characterized by the number of times that the baryonic 
fluid cycles between ingoing and outgoing. If the conductivity is at the continuum threshold k^, 
then the solution cycles repeatedly between ingoing and outgoing, displaying a discrete self-similarity 
reminiscent of that observed in critical collapse. Below the continuum threshold k™, and except 
at the discrete values k„, the baryonic fluid drops through the Cauchy horizon, and in this case 
undergoes a shock, downstream of which the solution terminates at an irregular sonic point where 
the proper acceleration diverges, and there is no consistent self-similar continuation to zero radius. 
As far as the solution can be followed inside the Cauchy horizon, the radial direction is timelike. 
If the radial direction remains timelike to zero radius (which cannot be confirmed because the self- 
similar solutions terminate), then there is presumably a spacelike singularity at zero radius inside 
the Cauchy horizon, which is distinctly different from the vacuum (Reissner-Nordstrom) solution 
for a charged black hole. 

PACS numbers: 04.20.-q 



I. INTRODUCTION 



What really happens inside black holes? Despite 
substantial progress, particularly in the last decade and 
a half since the landmark papers by Poisson & Israel 
(1990) j6|J on mass inflation at the Cauchy horizon, 
by Ori & Piran (1990) [64[ on similarity solutions in 
general relativistic collapse, and by Choptuik (1993) 
on critical collapse and discrete self-similarity, the answer 
to this question remains unresolved, e.g. |3l l2ril2^ . l3rill60| 
and references therein. 

It is well known that the vacuum solutions for 
charged (Reissner-Nordstrom) and rotating (Kerr- 
Newman) black holes are not physically consistent as 
endpoints of realistic gravitational collapse, because their 
cores are gravitationally repulsive. For a charged black 
hole, the gravitational repulsion comes from the negative 
radial pressure (radial tension) of the electric field, while 
for a rotating black hole, the gravitational repulsion 
comes from the centrifugal force. The gravitational 
repulsion causes accreted material to tend to pile up 
in the subluminal region inside the inner horizon of 
the vacuum black hole, rather than falling on to the 
singularity, contradicting the proposition that the black 
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hole is empty outside the singularity. Only in the 
unique case of an uncharged, non-rotating black hole 
(Schwarzschild) is the vacuum solution a consistent 
endpoint of realistic gravitational collapse. 

The present paper follows a well-traveled trail of 
exploration into the interiors of black holes, by consid- 
ering only spherically symmetric self-similar solutions, 
and imagining that it might be reasonable to regard 
charge as a surrogate for rotation, e.g. p^ . General 
relativistic spherically symmetric self-similar solutions 
were first considered by Cahill & Taub (1971) ^(|, and 
were first applied to the problem of general relativistic 
gravitational collapse by Ori and Piran (1990) [b3 |. 
following earlier work on general relativistic self-similar 
solutions in cosmology. Since then there have been many 
investigations; see especially the review s [l7t 4f3 and 
references therein; and more recently pH l42l 157 . l59l Ir37l 

MM 

The driving goal of this paper is not so much to 
study the formation of a black hole by gravitational 
collapse, but rather to explore the interior structure of 
black holes after their formation. We have in mind the 
situation of a realistic astronomical black hole, perhaps 
stellar-sized, perhaps supermassive, which is being fed 
by accretion of baryonic matter. The accretion rates 
in astronomical black holes are typically tiny, in the 
sense that the accretion timescale, which might be 
millions to billions of years, is vast compared to the 
characteristic timescale of the black hole, which might be 
milliseconds to hours. Thus it is to be expected that the 
superluminal region outside the inner horizon of a black 
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hole should closely approximate the vacuum solution. 
The interesting question is where and how the interior 
structure of the black hole deviates from the vacuum 
solution. 

In a seminal paper, Poisson & Israel (1990) showed 
that if ingoing and outgoing fluids are allowed to pass 
through each other inside a charged black hole, then the 
generic consequence is 'mass inflation' as the counter- 
streaming fluids approach an inner horizon. During 
mass inflatiom the interior mass, the so-called Misner- 
Sharp mass |58| . a gauge-invariant scalar quantity, 
exponentiates to enormous values. The phenomenon 
of mass inflation has been confirmed a naly tically and 
numerically in many papers 0, la, LLil LUJ UM 113 LiJ ■ 

Physically, mass inflation is a consequence of the fact 
that, if ingoing and outgoing fluids drop through an inner 
horizon, then they must necessarily pass into causally 
disconnected parts of spacetime. The only way this 
can happen is for the ingoing and outgoing fluids to 
exceed the speed of light relative to each other, which is 
impossible. As the counter-streaming fluids race through 
each other at ever closer to the speed of light in their 
attempt to drop through the inner horizon, they generate 
a large radial pressure. The growing radial pressure 
amplifies the gravitational force that accelerates the 
ingoing and outgoing fluids through each other, which 
results in mass inflation (see §V of [37j for a more 
comprehensive discussion). 

Analytic and numerical work on spherically symmetric 
collapse has commonly modeled the fluid accreted by the 
black hole as a massless scalar field, usually taken to be 
uncharged i IH E El ElJllMMllJllil M 
H|, but sometimes charged |26l HfSl Idl III. l6lt l7lj. A 
key property of a massless scalar field is that it allows 
waves to counter-stream relativistically through each 
other, which allows the phenomenon of mass inflation 
on the Cauchy horizon to occur. In the present paper 
and its companion, hereafter Paper 2 |33, we choose to 
adopt a somewhat different point of view, treating it as a 
matter of physics as to whether or not to allow relativistic 
counter-streaming, rather than assuming a fluid that has 
that property built in. 

In this paper, we assume that the black hole is fed in a 
self-similar fashion with a charged, conducting, baryonic 
fluid having a relativistic equation of state Pb/pb = 1/3. 
Although the fluid is electrically conducting, so that 
oppositely charged fluids do drift through each other, this 
proves not enough to permit mass inflation. The charged 
baryonic fluid is electrically repelled by the charge 
of the black hole, generated self-consistently by the 
charge of the accreted matter, and if the conductivity is 
small enough, then the baryonic fluid naturally becomes 
outgoing, and can drop through the outgoing inner 
horizon, the Cauchy horizon. 

The electrical conductivity of the charged baryonic 
fluid is treated in this paper as an adjustable parameter. 
If the conductivity is set at a realistic level, then the black 
hole accretes only if the charge-to-mass ratio is set to a 



tiny value. This is of course a consequence of the huge 
charge-to-mass ratio of individual protons and electrons, 
e/m p w 10 18 , where e is the dimensionless charge of the 
proton or electron, the square root of the fine-structure 
constant, and m p is the proton mass in units of the 
Planck mass. However, electrical conduction in a charged 
black hole is analogous to angular momentum transport 
in a rotating black hole, and angular momentum 
transport is a much weaker process than electrical 
conduction (if angular momentum transport were as 
strong as electrical conduction, then accretion disks 
would shed angular momentum as quickly as they shed 
charge, and accretion disks would not rotate). In 
the interests of regarding charge as a surrogate for 
angular momentum, it makes sense to investigate how 
the conductivity affects the internal structure of the black 
hole, with the conductivity being greatly suppressed 
compared to any realistic conductivity, but nevertheless 
possibly consistent with what might be a reasonable 
rate for the analogous angular momentum transport in a 
rotating black hole. 

In Paper 2 [37J we allow the black hole to accrete, 
in addition to the charged baryonic fluid, a pressureless 
neutral 'dark matter' fluid. The point of introducing dark 
matter is specifically to permit mass inflation, caused by 
relativistic counter-streaming between outgoing baryons 
and ingoing dark matter near the inner horizon. We defer 
further discussion of mass inflation to Paper 2. 

Previous investigations of self-similar solutions in 
gravitational collapse have commonly required that the 
solutions be reg ular at the origin, where the black hole 
first forms 1171 |3J, HJ ||. The present paper does 
not impose this requirement from the outset, but rather 
establishes boundary conditions outside the black hole, 
and integrates the equations inward. The outer boundary 
is set at the sonic point where the infalling baryonic 
fluid is assumed to transition smoothly from subsonic 
to supersonic velocity. This point of view seems natural, 
since information can only propagate inwards inside the 
black hole. 

A final difference between the present paper and earlier 
works is mathematical rather than physical. In this paper 
and its companion we choose to use the tetrad formalism, 
e -g- [Tol I22I l3lj . While it is to some extent a matter of 
taste what gauge or formalism to adopt, physical results 
being independent of such choices, the tetrad formalism 
does reveal the physics in an elegant and transparent 
way. We owe much to the opus by Lasenby et al. [53| . 
which showed us how to proceed at an early stage of this 
project. In particular, the gauge choices adopted here, 
equation JSJ, are those suggested by |53|. 

The structure of this paper is as follows. Section [n] 
presents the general relativistic equations governing the 
interior and exterior structure of a spherically symmetric 
black hole accreting a charged, electrically conducting 
fluid of baryons. Section IIIII introduces the hypothesis 
of self-similarity, and sets out the equations that follow 
from that hypothesis. Section llVl gives results for black 
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holes that accrete only baryons. Section Ivl addresses the 
question of what it would actually look like if you fell into 
one of the black holes described in this paper. Finally, 
section IVll summarizes the findings of this paper. 



II. EQUATIONS 

This section presents the general relativistic equations 
governing a spherically symmetric black hole accreting a 
charged, electrically conducting fluid of baryons. 



Tetrad formalism 



Let £ M denote a system of spacetime coordinates, 
and let denote the basis of tangent vectors in that 
coordinate system. By definition, the scalar products of 
the tangent vectors constitute the metric g^ v 



a^9u = g^ v 



(i) 



In the orthonormal tetrad formalism, a set of locally 
inertial frames is erected at each point of the spacetime. 
The locally inertial frame at each point has axes 7 m , the 
tetrad, whose dot products form, by construction, the 
Minkowski metric rj m n 



7m ' in Vn 



(2) 



In this paper dummy latin indices signify locally inertial, 
or tetrad, frames, while dummy greek indices signify 
curved coordinate frames. The axes 7 m of the locally 
inertial frames are related to the tangent vectors by 
the vierbein e„, M and its inverse e m „ 



7» 



9{i — 6 mTd 



(3) 



The vierbein provide the means of transforming between 
the tetrad components and coordinate components of 
any 4-vector or tensor object. For example, the 
tetrad components p m and the coordinate components 



p M of the 4-vector p 



P m lm are related 



by p m = e m M p M . Tetrad components are raised, 
lowered, and contracted with the Minkowski metric ?y mn , 
while coordinate components are raised, lowered, and 
contracted with the coordinate metric g^. 

The most general form of the vierbein e m M consistent 
with spherical symmetry is |69|, in Cartesian coordinates, 



eo 



a , e. 



eo' 



jXiXj + A (Sij - Xi&j) + fie ljk x k (4) 



where is the time index, i,j, k are taken to run 
over spatial indices 1,2,3, and Sijk is the completely 
antisymmetric flat space spatial 3-tensor. Appendix [S] 
gives results for the general case, but here it is convenient 
immediately to make the gauge choices 



, A = 1 , /i = 



(5) 



so that the vierbein e m M simplify to 

eo° = a , — , e l — (3x l , 

e-i = 'yXiXj + (S^ - XiXj) . (6) 

The gauge choices v = and /i = can be effected 
by starting from a general vierbein of the form 
and respectively boosting and rotating the tetrad frame 
appropriately at each point. The gauge choice A = 1 then 
corresponds to scaling the radial coordinate r so that 
the proper circumference of a circle of radius r is 2irr 1 a 
traditional and natural choice. The inverse vierbein e m ^ 
corresponding to the vierbein of equations © is 



e°o = i, e° 4 = 



e o 



(7) 



Directed derivatives d m are defined to be the spacetime 
derivatives along the axes j m of the tetrad frame: 



d 

f) = <v -Be v 



(8) 



where d = g^d/dx^ = g^g v d/dx^ is the invariant 
spacetime vector derivative. The directed derivatives 
d m depend only on the choice of tetrad frame, and are 
independent of the choice of coordinate system. Unlike 
the coordinate partial derivatives d/dx^ 1 , the directed 
derivatives d m do not commute. For the present paper, 
the most important directed derivatives are the directed 
time derivative dt and the directed radial derivative d r , 
which for the vierbein of equations © are 



dt 
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Tt 



d_ 

dr 



r 
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d_ 

dr 



(9) 



The directed time derivative dt is the proper time 
derivative (commonly written d/dr) along the worldlinc 
of an object instantaneously at rest in the tetrad frame. 
The coordinate 4-velocity u M = dtx^ of an object 
instantaneously at rest in the tetrad frame is given by 
the proper time derivatives of the coordinate time t and 
radius r 



dtt 



d t r = p 



(10) 



The directed radial derivative d r is the proper radial 
derivative from the point of view of an object instan- 
taneously at rest in the tetrad frame. The proper radial 
derivatives of the coordinate time t and radius r are 



d r t = , 8 r r = 7 



(11) 



For the more general vierbein of equation the 
directed radial derivative would be d r = vd/dt + ^d/dr, 
and it seen that the gauge choice v = is equivalent to 
choosing the time coordinate t so that the proper radial 
derivative d r is proportional to the coordinate radial 
derivative d/dr. 
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It is useful to note that the vierbein coefficients (/?, 7) 
constitute the time and radial components of a covariant 
tetrad frame 4-vector, the radial 4-gradient d m r 



(3 = d t r , 7 = d r i 



(12) 



It follows that /3 2 — 7 2 is a Lorentz scalar. Moreover, 
j3 2 — 7 2 is gauge-invariant with respect to arbitrary 
transformations of the time coordinate t — > t(t',r), the 
only remaining gauge freedom admitted by the vierbein 
of equations ©. The gauge-invariant scalar f3 2 — 7 2 is 
related to the mass M interior to r, equation Q24J) below. 

The vierbein coefficient a = dtt is not the time 
component of a tetrad frame 4-vector because, unlike the 
radial coordinate r, the time coordinate t changes when 
the tetrad frame is boosted, the gauge of time being fixed 
in accordance with equations In flIII Bl it will be seen 
that, for self-similar solutions, the quantity t/(ar) = £ 
is the time component of the homothetic Killing 4-vector 

The metric ds 2 = g^dx^dx" — r) mn e m p e" v dx^ dx v 
corresponding to the vierbein of equations JHJ is 



ds 2 = - 



d,t 2 



1 

~2 



dr 



a J 



+ r 2 do 2 



(13) 



where do 2 = d0 2 + sin 2 6 d(f> 2 is the metric of the surface 
of a unit 3-sphere. It is apparent that the vierbein 
coefficients a and (5 are related to the lapse and shift 
in the ADM formalism Q (see e.g. [&4| for a pedagogical 
review), the lapse being 1/a, and the shift being —(3/a. 
The reciprocal coefficient 1 /j might be termed the radial 
stretch, since the proper radial interval at fixed time is 
dr/j. 

An object instantaneously at rest in the tetrad frame 
satisfies dr/dt = j3/a, according to equation l|10|l . and it 
then follows from the metric that the proper time 
experienced by the object is dt/a. 

In the tetrad formalism, the covariant derivative D n of 
the tetrad components p k of a 4-vector is 



D n k — 3 n k + T fe »" 



(14) 



where the tetrad frame connection coefficients T k nnl also 
known as the Ricci rotation coefficients, are defined, at 
least in the case of vanishing torsion, by the directed 
derivatives of the tetrad axes 



d nlm = T« mnlk . (15) 
In terms of the vierbein derivatives dkmn defined by 



Bp a 



(16) 



the connection coefficients Tkmn = Vkl^mn w hh all 
indices lowered are 



km a 
1 



(17) 



(dk 



dmkn ~t~ d nra k d n fc m -)- d mn fc dknra) 



The tetrad frame connection coefficients T^nm are 
antisymmetric in their first two indices, 

Tfcmn — f mfcn (-^^) 

which expresses mathematically the fact that T k nn for 
each given final index n is the generator of a Lorentz 
transformation between tetrad frames. The usual 
coordinate frame connection coefficients, the Christoffel 
symbols T Klll , = <? K Ar^„, are related to the tetrad frame 
connection coefficients T^nm by 



T — e k p m e n (Yi - d 



kmn J 



(19) 



which has the usual T 



n^u — r K y M symmetry for vanishing 
torsion. 

For the vierbein given by equations JSJ, the non-zero 
tetrad frame connection coefficients are, in Cartesian 
coordinates [coefficients with switched first two indices 
follow from antisymmetry, eq. (|18|) ]. 



iOO 



- iOj 



ijk 



gxi 



7-1 



(20a) 
(20b) 

(20c) 



where g in equation l|20a|) is the proper acceleration 
experienced by an object that is comoving with the tetrad 
frame 



d r In a 



(21) 



and h in equation ij20bfl is the proper radial gradient 
of the proper velocity between objects each of which is 
comoving with the tetrad frame at its position 



d(3 dlna 

h= dthij-P — — 

or or 



(22) 



The Riemann tensor Rkimn in the tetrad frame is 
defined in the usual way by the commutator of the 
covariant derivative, [D^, Di]p m = RkimnP n , and is given 
in terms of the tetrad frame connection coefficients by 

Rklmn — &l^nmk &k^nrnl ^mk^ n 3^ ^ml^ n jk 

+ (r^j, — T J kl )r nm j (23) 

which has two extra terms (the last two) compared to 
the usual coordinate expression for the Riemann tensor 
in terms of Christoffel symbols. The Ricci tensor and 
scalar are given by the usual contractions of the Riemann 
tensor, R km = T] ln R klmn and R = r/ km R km , and the 
Einstein tensor is given by the usual expression in terms 
of the Ricci tensor and scalar, Gk m = Rkm — \Rf]km- 

In the present case where the vierbein are given by 
equation ijBJ), define the interior mass M, the Misner- 
Sharp mass, by 



2M 
r 



1 + (3 2 - 7 2 



(24) 
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As remarked following equations (|12f> . 1 — 7 2 is a scalar, 
a gauge- invariant quantity, and thus so also is the interior 
mass M. Further, define the symbols F, R (not the Ricci 
scalar!), P, and S by 



F 
R 
P 
2S 



l(da-0g) = ^(drP-jh) , (25a) 



f ( 1 



7 \ r 



T d r M-(3F) , 



r \ J 
d r P + g(R + P) 



(25b) 
(25c) 







■d t F + 2[ ^ + h)F 
r 



In terms of these quantities, the components of the 
Einstein tensor G mn in the tetrad frame are, in Cartesian 
coordinates, 

G 00 = 2R , (26a) 
G m = -2F Xi , (26b) 
Gij = 2PS lJ +2S(S lJ -x l x J ) . (26c) 

The Einstein tensor embodies the compressive part of the 
Riemann tensor. The non-compressive, or tidal, part of 
the Riemann tensor is the Weyl tensor Ckimn , which here 
simplifies to 

Ckimn = G [tkltmn %kl&mn g {l]km^]ln ^7/cn^/m)] 

(27) 



where C is 



C = R + S 



3M 



(28) 



and the bivectors (antisymmetric tensors) t mn and x Ti 
have non-zero components 



iO 



-toi = Xi 



^ijk^k 



(29) 



B. Einstein equations 



The Einstein equations in the tetrad frame are 



G r . 



8ttT„ 



(30) 



where T mn is the energy-momentum tensor. It is 
apparent from equations l|26(l for the Einstein tensor G mn 
that the Einstein equations will take their simplest form 
when expressed in the center of mass frame, defined to 
be the frame in which the momentum density vanishes, 
= 0. In the center of mass frame, Gio vanishes, and 
equation (|26bfl then implies that the quantity F defined 
by equation (|25a|l vanishes, 



F = 



(31) 



significantly simplifying the expressions for the other 
components of the Einstein tensor. In the center of mass 
frame, the most general form of the energy-momentum 
tensor T mn consistent with spherical symmetry is 



Ti, 



p r x l x j + p± (<% - 



(32) 



where p is the proper energy density, p r the proper radial 
pressure, and p± the proper transverse pressure. The 
expressions (|26ll for the Einstein tensor and (|32() for 
the energy-momentum tensor inserted into the Einstein 
equations (|30|l then imply 



R = Anp , P = Anp r , S = 47r(p ± - p r ) 



(25d) The resulting Einstein equations are 



dM 
dr 

d tl - p g = o , 



Airr 2 p 



dtP -19 + ~2 
27, 



+ 4irrp r = , 



d r Pr (j)± - Pr) + {p + Pr)9 = 

r 



(33) 

(34a) 
(34b) 
(34c) 

(34d) 



Equations l|34l) , coupled with the defining equations 121|) 
for the acceleration g and (|24(l for the interior mass 
M, constitute the Einstein equations in the most 
general form consistent with spherical symmetry. The 
first and last equations may be regarded as elliptic 
constraint equations, the first equation 134afl expressing 
what appears to be the familiar relation between mass 
and density, and the last equation 134dl) expressing 
pressure balance, the Euler equation. The midde two 
equations are hyperbolic evolution equations governing 
the evolution of the metric coefficients 7 and f3. 
Equation (|34c|l relates the acceleration d t f3 of the tetrad 
frame to the gravitational force, which comprises the 
familiar Newtonian force M/r 2 , minus a term 717 that 
comes from the acceleration generated as a result of 
pressure balance, plus a term Airrp r proportional to the 
radial pressure p r . 

It is this last term Airrp r in the acceleration equa- 
tion (I34c|) that makes the interiors of charged black holes 
problematic and interesting, for it is this term that, 
thanks to the negative radial pressure of the electric 
field, tends to make charged black holes gravitationally 
repulsive in their deep interiors. 

The result of taking (3 times equation l|34cf) minus 7 
times equation l|34b|) is 



(35) 



d t M + 4Trr 2 p r (3 = 



which expresses energy conservation, the first law of 
thermodynamics (applied to the baryonic fluid, not to the 
black hole as a whole). Equations (|34a|l and (|35|l amply 
justify the identification of M defined by equation l|24|) 
as the mass (energy) interior to radius r. 
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Conservation of energy-momentum as expressed by 
the vanishing of the covariant derivative of the energy- 
momentum tensor, 



D m T mn = 



(36) 



is automatically built into Einstein's equations, a 
consequence of the Bianchi identities, which ensure that 
D m G mn = 0. The time component, n = 0, of 
equation l|3*6"jl gives the energy conservation equation 



Otp-\ {p± 

r 



Pr) + 



(p + p r )dr 2 f3 
r 2 dr 







(37) 



which can also be derived from Einstein's equations l|34(l 
by taking the time derivative dt of equation (|34a|) . elim- 
inating dtM using equation (|35|l . and then simplifying 
using the Euler equation i|34d|) . The spatial components, 
n = 1,2,3, of equation H3fi|l give a momentum 
conservation equation which reduces precisely to the 
Euler equation (|34d(l . 

The familiar Oppcnheimcr-Volkov equation for general 
relativistic hydrostatic equilibrium 



dp (p + p) (M + 47rr 3 p) 
= r 2 (l - 2M/r) 



(38) 



is recovered from the Einstein equations Q34[) by setting 
the tetrad frame at rest, equivalent to setting (3 = dtX — 
0, by further assuming isotropic pressure, p r = p±, and 
finally by eliminating the acceleration g in the Euler 
equation i)34dl) using the acceleration equation (|34c|l . 



C. Baryonic fluid 

As noted by [H, [H HE El and many others, the 
Einstein equations admit similarity solutions only if 
pressure and density are proportional. In this paper we 
assume that the black hole accretes a baryonic fluid with 
an isotropic pressure and a relativistic equation of state 



Pb = wp b 



w 



(39) 



In reality, the baryonic fluid accreting on to astronomical 
black holes is likely to be near but not fully relativistic. 
Nevertheless, given that self-similarity requires w to 
be constant, the choice w = 1/3 seems the most 
reasonable. Moreover, it may be anticipated that the 
fluid will be adiabatically compressed by gravitational 
repulsion inside the black hole, and that it may even 
undergo a shock if it passes through an inner horizon. 
Both processes will render the baryonic fluid yet more 
relativistic. 



D. Electromagnetic field 

The assumption of spherical symmetry implies that 
the electromagnetic field can consist only of a radial 



electric field Q/r 2 . The only non-zero components of 
the electromagnetic field tensor F mn in the tetrad frame 
are then 



F, 



i0 



-Foi — —^Xi . 



(40) 



Equation (|40|l may be regarded as defining what is meant 
by the charge Q interior to radius r. The electromagnetic 
energy-momentum tensor in the tetrad frame is given by 
4nT mn = F mk F n k - \ri mn F k iF kl , which implies that the 
density p e and radial and transverse pressures p e>r and 
p Cy ± of the electric field are 



-Pe 



: Pe 



8-Kr 4 



(41) 



Notice that the radial pressure p etT is negative; it is 
this radial tension of the electric field that leads to 
gravitational repulsion inside charged black holes. 

An electromagnetic field that consists only of a radial 
electric field automatically satisfies two of Maxwell's 
equations, the source-free ones. The other two Maxwell's 
equations, the source ones, are, in the tetrad frame, 



D„F mn = d n F n 



1 nl^ 



+ r; d F ml = A7rf 



(42) 



For the electromagnetic field given by equation 14Q|I , and 
with j° = q, j % — the Maxwell equations l|42|) reduce 
to 



d r Q — 47rr q , 
D t Q = -ATtr 2 j 



(43a) 
(43b) 



The quantity q is the proper charge density in the fluid 
frame, while j is the proper radial current. 

There is no need to adjoin the Lorentz force law, 
because that is built into Einstein's equations, which 
automatically enforce energy-momentum conservation. 
Specifically, if the density and pressure in the Euler 
equation l|34d|) are written as a sum of baryonic 
and electromagnetic contributions, the electromagnetic 
pressure and density being given by equations l|41|l . 
and if the resulting radial derivative of charge Q is 
eliminated using the Maxwell equation (|43a|l . then the 
Euler equation (|34d|) for the charged baryonic fluid 
becomes 



(44) 



and it is apparent that the qQ/r 2 term expresses the 
Lorentz force. 

Similarly, electrical conduction causes Ohmic genera- 
tion of heat, but there is no need to doctor the Einstein 
equations since they already automatically include this 
effect. Specifically, the electromagnetic contribution to 
dtp in the energy conservation equation (|37|l involves a 
term ~ (QdtQ) /r 4 ~ jQ/r 2 which is precisely the Ohmic 
dissipation term. 
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E. Conductivity 

It is still necessary to specify a governing equation for 
the radial current j. For diffusive electrical conduction, 
the radial current j is proportional to the radial electric 
field Q jr 2 , with the coefficient of proportionality defining 
the electrical conductivity a 



Q 

3 = °7> 



(45) 



which is just Ohm's law. A realistic value (which is not 
used in this paper) of the electrical conductivity of a 
baryonic plasma at a relativistic temperature T is [l| 



C kT 
e 2 lne _1 h 



(46) 



where e is the dimensionless charge of the electron, 
the square root of the fine-structure constant, and the 
factor C » 15 depends on the mix of particle species. 
This electrical conductivity is huge. A dimensionless 
measure of the conductivity (which has units 1/time) 
is the conductivity a times the characteristic timescale 
= GM/c 3 of the black hole, which is of order 



at 



BH 



T 



(47) 



where fcXeH = S/^bh is the characteristic temperature 
of the black hole (for a Schwarzschild black hole, this 
characteristic temperature Tbh is 87r times the Hawking 
temperature). In the astronomical situation considered 
here the temperature T of the plasma is huge compared 
to the characteristic temperature Tbh of the black hole. 
Indeed if this were not so, then mass loss by Hawking 
radiation would tend to compete with mass gain by 
accretion, an entirely different situation from the one 
envisaged here. 

In the present paper, we do not use a realistic value 
for the conductivity, equation (|46[) . but instead adopt a 
phenomenological conductivity which is greatly reduced 
compared to any realistic value. As remarked in the 
introduction, the point is that if charge is regarded 
as a surrogate for angular momentum, then electrical 
conduction can be considered as an analog of angular 
momentum transport, which is intrinsically a much 
weaker process than electrical conduction. If it is 
assumed that the phenomenological conductivity a is 
some function of the baryonic density pb, then the 
hypothesis of self-similarity requires that this function 
be the square root (dimensional analysis shows that ar 
is dimensionless, and p^r 2 is dimensionless) so 



1/2 

(4tt) 1 /2 



(48) 



where k is a phenomenological dimensionless coefficient 
of conductivity, the factor of l/(47r) 1 / 2 in equation l|48|) 
being inserted to simplify the equation when cast in 



self-similar form, equation (|51|l below. The realistic 
conductivity given by equation (|46|l is approximately 

proportional to p^ 4 , so the notion that the conductivity 
increases as some positive power of baryonic density 
seems physically reasonable. 

It should be noted that the radial current is expected 
to increase slightly the radial pressure compared to the 
transverse pressure of the charged baryonic fluid, but 
this increase is small for diffusive electrical conduction 
(and utterly miniscule with the suppressed conductivity 
adopted here), and is neglected in this paper. 



III. SIMILARITY SOLUTIONS 

A. Similarity hypothesis 

Dimensional analysis of the spherically symmetric 
coupled Einstein-Maxwell equations l|34|) and i|43[l . 
along with the definitions i|21|) of the acceleration 
g and l|24|) of the interior mass M, combined with 
the equation of state (|39Jl of the baryonic fluid, the 
energy-momentum l|41l) of the electric field, and the 
phenomenological conductivity l|48l) , reveals the following 
quantities to be dimensionless: 



_ or 

y = gr , z = Airr" pb , 



7, ^, 2, (49 ) 
r r 

z„ = 47rr 2 q , s = Aura . 



It is convenient also to denote the dimensionless energy 
density 47rr 2 p e of the electric field by 



A 1 

z e = 4^ p e = — . 



(50) 



For the phenomenological conductivity a given by 
equation Q48[l. the dimensionless conductivity s is 



a = nz 1 ' 2 



(51) 



If further the gauge of the time coordinate t is fixed 
in a natural way, by setting t equal to the proper time 
experienced by an object at some boundary, then the 
dimension of time t is the same as the dimension of radius 
r, and then a and r/t are separately dimensionless, not 
just their product 77 = ar/t. In this paper it is not really 
necessary to choose a specific time coordinate t, because 
no observable quantities depend on the choice. All the 
same, there is a natural choice available if needed, which 
is to synchronize the time to the proper time recorded by 
objects (dark matter test particles) that free-fall radially 
starting from zero velocity far away. 

The dimensionless variables of equation (|49|) form a 
(more than) complete set for the problem at hand, and 
it follows [IT} that the spherically symmetric Einstein- 
Maxwell and subsidiary equations admit similarity 
solutions in which the dimensionless variables are all 
functions of the dimensionless variable r/t. In particular, 
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the mass M, the charge Q, and the radius r of the black 
hole, evaluated at some similarity point, such as the outer 
horizon, increase linearly with time 



M oc Q oc r oc t 



(52) 



Actually the Einstein-Maxwell and subsidiary equa- 
tions would seem to admit more general solutions in 
which the dimcnsionless variables given by equation 149(1 
are considered to be functions of r/F(t) with F(t) 
some arbitrary function of time t. However, this 
apparently greater generality is not genuine, but merely 
an artifact of the gauge freedom in the choice of the 
time coordinate t. Nevertheless, the gauge freedom does 
have the consquence that the similarity equations do not 
involve a and r/t separately, but only their product, the 
dimcnsionless variable rj = ar/t. 

Spacetime points that are hxed in the self-similar frame 
have coordinate velocity v r /v t — r/t — constant, where 
f M is the 4-velocity in the coordinate frame. In the 
baryonic tetrad frame these points move with radial 
velocity V = u r /u l , where u m = e m ^v^ is the 4-velocity 
in the baryonic tetrad frame (the 4-velocity u m in the 
tetrad frames is denoted by a latin u, whereas the 4- 
velocity in the coordinate frame is denoted by a 
greek upsilon). With the inverse vierbein e m M given by 
equations (7J , the radial velocity V of the similarity frame 
relative to the baryonic tetrad frame is 



V = 



v-P 
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(53) 



This velocity V is minus the proper velocity of the 
baryonic fluid through the self-similar frame. Horizons 
occur where the velocity equals the speed of light 



V = ±1 at horizons 



(54) 



while sonic points occur where the velocity equals the 
sound velocity 



V = ±^/w at sonic points 



B. Homothetic Killing vector 



(55) 



Similarity solutions are invariant under a scale trans- 
formation of t and r at hxed r/t. The generator of this 
scale transformation is the homothetic Killing vector £ 
given by H E E3 IH |H E H H 



. _ d d d 

rE • O — t — r — — 

dt dr dhit 



(56) 



r/t 



The three right hand sides of equation (|56[1 express the 
homothetic vector £ in three frames: the (t, r) coordinate 
frame; another coordinate frame in which the time and 
radial coordinates are ln£ and r/t; and the baryonic 
tetrad frame. The components of the homothetic 4- 
vector £ in the (t,r) coordinate frame are x^/r, with 



time and radial components (t/r, 1). Equation l|56(l 
coupled with the relations © between directed and 
coordinate derivatives imply that the components £ m of 
the homothetic 4- vector in the baryonic tetrad frame are 



e = 



V 
V 



(57) 



where 77 = ar/t, and the velocity V is given by 
equation l|53l) . In terms of the homothetic vector, the 
relation (|53[1 for the velocity V can be re-expressed as 
the statement that the scalar product of the covariant 4- 
vector (/3, 7) with the contravariant 4-vector £ m equals 
one, + 7^'" = 1. The magnitude squared of the 
homothetic vector defines the homothetic scalar H 



H 



1-V 2 
V 2 



(58) 



which, like the interior mass M, is gauge- invariant. 
Horizons occur where the homothetic scalar vanishes, 
H = 0, which demonstrates that the location of horizons 
is independent of the choice of the coordinate system or 
of the frame of reference, which is as it should be. 

Associated with the scale transformation symmetry 
is a conservation law. To derive this conservation law, 
consider the Lagrangian for a (neutral) test particle 
freely- falling with coordinate 4-velocity v^ 1 = dx^/dr 



£ = 



(59) 



If the metric (|13fl is cast in terms of coordinates lni 
and r/t in place of t and r, then all of the metric 
coefficients become proportional to t 2 at fixed r/t, so 
that dC/d\nt — 2C — —2, the last equality being true at 
least for a massive particle, for which v^v^ = — 1. The 
equation of motion for the coordinate In t of the freely- 
falling particle is 



d dC 
dr dv lnt ~ 

which then reduces to 

dvint 
dr 

Equation (|61|l integrates to 
Vint + 



dC 
dhit 



= 



1 = 







(60) 



(61) 



(62) 



where the constant of integration has been absorbed 
into a shift of the zero point of the proper time r. 
Equation (|62|l is the sought-for conservation law: it says 
that the homothetic momentum v\ nt associated with the 
coordinate hit of a freely-falling particle is equal to minus 
the proper time r along the geodesic, provided that the 
zero point of proper time is set appropriately. In the 
tetrad frame, the homothetic momentum v\ n t is 



Vint = r £ m u r ' 



(63) 
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Although the above derivation was for a massive particle, 
the results <|62[) and (|63l) remain valid also for a massless 
particle, with the understanding that for a massless 
particle the proper time r is a constant along the null 
geodesic. 



C. Ingoing vs. Outgoing 

In the Reissner-Nordstrom (RN) geometry, geodesies 
between the outer and inner horizons can be classified as 
ingoing or outgoing according to the sign of their specific 
energy E = —Vt, where w M is the coordinate 4- velocity 
along the geodesic. Outgoing geodesies have negative 
energy, and go backwards in RN time t. The horizons 
of the RN geometry can likewise be classified as ingoing 
or outgoing. Only ingoing geodesies (those with E > 0) 
can pass through ingoing horizons, while only outgoing 
geodesies (those with E < 0) can pass through outgoing 
horizons. Timelike or null geodesies that start from 
outside the outer horizon necessarily go forward in time 
and have positive energy, and are therefore necessarily 
ingoing, and cannot pass through the outgoing inner 
horizon, the so-called Cauchy horizon. Charged particles 
do not follow geodesies, being accelerated by the electric 
field, and particles with the same sign of charge as the 
black hole, and with a large enough charge-to-mass ratio, 
do routinely switch from ingoing to outgoing between the 
outer and inner horizons, and thus do fall through the 
outgoing Cauchy horizon. 

In the similarity solutions, geodesies between the outer 
and inner horizons can be classified as ingoing or outgoing 
according to the sign of minus the homothetic momentum 
— Vint, equation (|63() . This definition of ingoing versus 
outgoing is gauge- invariant, independent of the choice of 
either coordinates or tetrad frame. Note that according 
to equation v\ nt can only decrease along a geodesic, 
since the proper time r can only increase. Thus a 
freely-falling (neutral) massive particle can change from 
outgoing (ui nt > 0) to ingoing (vi nt < 0), but once it is 
ingoing it must remain ingoing thereafter, as long as it 
continues in free-fall. The asymmetry between ingoing 
and outgoing reflects the fact that the similarity solution 
is not invariant with respect to time reversal. 

The tetrad frame can be classified as ingoing or 
outgoing according to whether objects instantaneously 
at rest in that frame are ingoing or outgoing. The 
4-velocity of an object instantaneously at rest in the 
tetrad frame is it* = 1 and u r — 0, and in this case 
equation 1|63|) simplifies to v\ n t = —r/rj. Thus the tetrad 
frame is ingoing or outgoing according to the sign of 
rj = ar/t. When the tetrad frame switches between 
ingoing and outgoing, so that rj changes sign, the time 
component £* = l/rj of the homothetic vector must pass 
through zero; it cannot pass through infinity because the 
homothetic vector £ m is finite. At the same time, the 
radial component C = V/rj must remain finite, neither 
zero nor infinite, again because the homothetic vector C n 



is finite, neither identically zero nor infinite. It follows 
that when the tetrad frame switches between ingoing and 
outgoing, rj must pass through ±oo, and the velocity V 
must simultaneously pass through ±oo, changing sign 
simultaneously with rj. Consistency of the self-similar 
solutions requires that C always remain of the same sign, 
since if C changed sign it would indicate that the fluid 
turns back on itself, which physically cannot happen (at 
least for baryons). In the black hole solutions of the 
present paper, C is initially positive (77 and V are both 
positive). Thus one concludes that the baryonic tetrad 
frame is ingoing if V is positive, and outgoing if V is 
negative, 



V > 

V < 



ingoing 
outgoing 



(64) 



and that the passage between ingoing and outgoing is 
marked by V passing through ±00. It is this condition on 
the velocity V that will be used in the results section HVI 
to indicate whether the baryonic fluid is ingoing or 
outgoing. 

There is no physical divergence or discontinuity in 
the properties of the baryonic fluid when the velocity 
V passes through ±00. All that happens is that the 
homothetic 4-vector £ in the baryonic tetrad frame 
switches between pointing forwards in time (ingoing, 
positive £') to backwards in time (outgoing, negative £*). 



D. Geodesies 

Towards the end of this paper, <|V] the visceral question 
is posed: What does it actually look like if you fall inside 
one of the black holes described in this paper? To address 
this question requires ray-tracing, which requires solving 
for geodesies, notably null geodesies. 

The equations of motion for freely falling test particles, 
massive or massless, take their most transparent form 
when expressed in self-similar, or homothetic, coordi- 
nates, where the homothetic symmetry is explicit. Let th 
denote the homothetic time, the time coordinate defined 
by the gauge choice JSJ when the tetrad frame is chosen to 
be comoving, possibly superluminally, with the similarity 
frame. Define the self-similar coordinate X by 



dx^ HdHr,th) 

7 e + pc 



(65) 



where H is the homothetic scalar, given by equation 15811 . 
The factor of H + PC) is included in equation ll6"5|) 
so that v = dX/dr, see equation (|68b|) below, remains 
finite and well-behaved everywhere along a geodesic. 
Note that 7^* + PC is a scalar, the scalar product of 
the covariant 4-vector (7, 0) with the contravariant 4- 
vector (C,C)- The metric with respect to homothetic 
coordinates In th and X is 



ds 2 = r 2 [ -Hdhxtl 



dX 2 



do' 



(66) 
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Without loss of generality, let a test particle move in 
the 9 = ir/2 plane, with varying azimuthal angle (f>. In 
homothetic coordinates In th and X, the 4- velocity 
satisfies 



Vlnt h 



,f = 



(67a) 
(67b) 
(67c) 



The first of equations l|67[) is equation (|62|l for the 
homothetic momentum, as derived in 3III Bl The second 
equation expresses conservation of angular momentum 
per unit energy J. The last equation expresses 
conservation of rest mass per unit energy n, which can 
be taken to be \i = 1 for a massive particle, or fi = for 
a massless particle. 

From equations l|67|l and the metric <|66|l it follows that 
in homothetic coordinates the coordinate 4-velocity 
of a particle with proper time r, angular momentum per 
unit energy J, and mass per unit energy [i is given by 



,lnt h 



,x 



Hr 2 



±-i[T 2 -tf(J 2 + A 2 )] 



J 



(68a) 
(68b) 
(68c) 



in which the sign of v x is determined according to 
whether the particle is moving outwards or inwards 
relative to the similarity frame. Relative to an observer 
at rest in the tetrad frame, the tetrad 4-velocity u m of 
the particle is 



t _ er±^ r [T 2 -H{J 2 +fi 2 r 2 )] 



1/2 



Hr 



(69a) 



I 1/2 



?r±? [r 2 -H{J 2 + ^ 2 r 2 )\ 
u = - — d (69b) 



Hr 



J 



(69c) 



where £ m are the components of the homothetic 4- vector 
in the tetrad frame, equation (|57l) . and the ± signs in the 
expressions for u* and u r are the same as the sign of v , 
equation JgSbJl. 

The shape of a geodesic trajectory can be obtained by 
integrating d<j)/dX — v^/v x . For photons, the proper 
time r is constant, and the rest mass \i is zero, and the 
equation for the angle <j> along a null geodesic reduces to 
an integral 



JdX 



(r 2 - HJ 2 ) 



1/2 



(70) 



A null geodesic passes through peri- or apo-apsis in 
the self-similar frame where the denominator of the 
integrand of equation l|70[) vanishes. The separatrix 
between null geodesies that do or do not fall into the 



black hole, the equivalent of the photon sphere, occurs 
where the denominator not only vanishes, but is a double 
root, which happens when (the subscript ph signifies the 
photon sphere equivalent) 



T 2 



dH 
dX 



= . 



(71) 



H=H„ 



An observer sees a particle, either a photon or a 
massive particle, to come from angle \ away from the 
direction to the center of the black hole. The angle 
X G [0, 7r] is given by 



tanx = — = 
u r 



j[?TTe [t 2 - H(J 2 + M V)] 1/2 } 



(£') 2 (J 2 +mV) 



, (72) 

in which the =F sign is opposite to the ± sign in the 
expression (|69bp for u r . 

The apparent edge of the black hole on the sky is set 
by photons from the photon sphere equivalent, whose 
proper time r p h and angular momentum J p h satisfy 
equations (|7T|l . The angular size x P h S [0,7r] of the black 
hole on the sky is then 



tan x P h 



(e) 2 



(73) 



in which the =p sign is negative for observers outside 
the radius of the photon sphere equivalent, r > r p h 
(the photons move outward in the similarity frame), and 
positive for observers inside, r < r p h (the photons move 
inward in the similarity frame). 

The observed blueshift of photons at the edge of the 
black hole equals the ratio w ph /u ph0 of the observed- 
to-emitted energy of a photon from the photon sphere 
equivalent. We take the emitted energy u ph to be the 
energy of the photon from the point of view of an observer 
at rest in the self-similar frame at the radius of the photon 
sphere equivalent, r p h. The observed blueshift factor is 
then 



» P h 

ph,0 



?"ph [H ph + (f 



r\21 



e<Tf (H ph -H) 



1/2 



(74) 



where again the =F sign is negative for observers outside 
the radius of the photon sphere equivalent, r > r p h, and 
positive for observers inside, r < r p h. 

In the limit of small accretion rates and small 
conductivities (and in the absence of mass inflation), 
the self-similar geometry asymptotes to the Reissner- 
Nordstrom geometry over any limited range of time, and 
X and H tend to 



X 



r z \ r 



(75) 



(76) 
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where 1 - 2M/r = 1 - 2M c /r + Q 2 /r 2 with M c and 
Q respectively the constant mass and charge of the RN 
black hole. In this limit, H(X) is a cubic polynomial 
in X for an uncharged black hole, Q = 0, or a quartic 
polynomial in X for a charged black hole, and the ray- 
tracing integral (j7t)(l becomes an elliptic integral, just as 
in the Schwarzschild and Reissner-Nordstrom geometries. 



E. Shock jump conditions 

The baryonic fluid may undergo a relativistic shock 
[H, Il6| . where the density, pressure, and velocity of the 
baryons change discontinuously across a shock front, a 
3-dimensional hypersurface in spacetime. Let n m denote 
the (spacelike) normal to the shock front. Conservation 
of energy-momentum imposes the shock jump conditions 



[n m T ri 



= 



(77) 



where []_ denotes the difference between post-shock 
(+) and pre-shock (— ) values. For a spherical shock 
wave, the shock normal n m in the rest frame of the 
shock front has components ■"' 
energy-momentum tensor T ( 
remains invariant under any radial Lorentz boost, so for 
a spherical shock wave the shock jump conditions reduce 
to jump conditions on the baryonic energy-momentum 



n- = 0, n r = 1. The 
of the radial electric field 



tensor T h mn alone 



r tt 



T l }L = . 



(78) 



Let V- and V+ denote the proper velocities of the 
shock front relative respectively to the pre-shock (— ) 
and post-shock (+) baryons. If the shock is self-similar, 
as considered in the present paper, then the velocities 
V± are equal to the velocities of the self-similar frame 
relative to the baryons, but equations (|79|) - (|82[) below 
remain valid also for a general spherical shock. The shock 
jump conditions (|78|l imply that the pre- and post-shock 
velocities V± are related by 

V-V+=w (79) 

and that the ratio of post- to pre-shock densities is 

Pb,+ V-(l 



Pb,- v+(\-vl) 



(80) 



Relative to the pre-shock baryons, the post-shock 
baryons are Lorentz boosted by 4-velocity u m given by 



(1-V-)(1 + V+) 



"I 1/2 



(1 + V-)(1-V+) 

1 



(81) 



equation 1|12|) . are Lorentz boosted by the 4-velocity 
u m across the shock. Written in a form that remains 
numerically well-behaved even under extreme conditions, 
the relation between post- and pre-shock coefficients is 



7 + ±/? + = (7-±/?-)(^±"'') 



(82) 



The homothetic Killing vector £ m , a contraviant 4- vector, 
is similarly Lorentz boosted across the shock: 



(£±C)(«*T« r ) 



(83) 



The scalars 1 — j 2 and H = — £ m £™ 1 are unchanged 
by a Lorentz boost, and are therefore continuous across 
the shock. The interior mass M, which is related to 
the scalar f3 2 — 7 2 by equation if^l , is consequently also 
continuous across the shock, as one might have expected. 
The interior charge Q is continuous across the shock, 
because the energy-momentum tensor of the electric field, 
equation 141|) . is invariant under a radial Lorentz boost. 



F. Integrals of the similarity equations 

The similarity hypothesis requires that all dimen- 
sionless quantities must be some function of a single 
dimensionless variable, which can be taken to be for 
example r/t. The directed time and radial derivatives 
of a dimensionless quantity f(r/t) are 



rd t f = (-r, + (3) 



df 



d\n(r/t) 



, rd r f = 7 



df 



d\n(r/t) 



(84) 



In particular, any dimensionless quantity f(r/t) must 
satisfy 



(d t + Vd r ) / = 



(85) 



where V, given by equation l|53|) . is the velocity of the 
self-similar frame relative to the tetrad frame. In terms of 
the homothetic Killing vector £ m given by equation 157fl , 
equation becomes 



rd m f = 



(86) 



The ordinary differential equations determining the 
self-similar evolution of the baryonic fluid admit three 
integrals. The first integral follows from 



rd m (M/r) = 



(87) 



as a particular case of equation l|8t)|) . Inserting equa- 
tions l|34a|l and l|35|l into equation l|87|) . and simplifying, 
yields an equation for the dimensionless ratio M/r of 
interior mass to radius 



M 



ziaC - «>/?£*) + 



(88) 



The vierbein coefficients (3 and 7, which form the 
time and space components of a covariant 4-vector, 



Equation (|88|l expresses M/r, defined in terms of /3 and 
7 by equation (|24|) . in terms of other dimensionsionless 
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variables. In the present paper we do not impose 
equation (|88|l as one of the evolution equations, but 
rather use it as a check on the accuracy of the numerical 
integration. Usually the check is satisfied to ~ 10~ 9 or 
better, but in some cases the check fails badly, namely 
in some cases where the integration terminates at an 
irregular sonic point below the Cauchy horizon. Where 
the check fails, the problem is that 7£ r — on the 

right hand side of equation l|88[l is a tiny difference of two 
large numbers, which leads to loss of precision. However, 
the problem is clearly with the check equation (|88|l . not 
with the equations being integrated, which is the reason 
for not including equation (|88|l among the evolution 
equations. 

The second integral of the similarity equations follows 
from 



rd m (Q/r) = 



(89) 



which, when the Maxwell's equations H43[) are inserted, 
yields an equation for the dimcnsionlcss charge density 



z q = 4irr 2 q 



0(1 + af) 



(90) 



The third integral of the similarity equations follows 
from 

w4Trr 2 (p + Pr )rd m l3 + ~fVrd m (4irr 2 Pr ) = (91) 

as another consequence of equation (|86|l . Equation l|91|l 
simplifies as follows. First recast the dtp r part of 
£ m d m p r in terms of dtPb and dtQ using the baryonic 
and electromagnetic equations of state and (|4"T)) . 

Then eliminate the derivatives dtfl, d r p r , dtp, and dtQ 
using equations l|34c[l . (|34d|l . (|3"T|) . and l|43b|) . along 
with the conductivity equation (|45|l . The result is an 
expression that contains no derivatives. Translating 
into the dimensionless variables of equations (|49|l and 
(|50|l yields an equation for the dimensionless proper 
acceleration y = gr 



ti 2wM/r + 2z e [(l-w) + (l-Ho)s£*] 

(i+ w )z[(e) 2 -^ 4 ) 2 ] 



(92) 



The denominator of expression l|92|l for the dimensionless 
acceleration y is proportional to V 2 — w, which is zero 
when the baryonic fluid velocity through the similarity 
frame equals the speed of sound, V — ±^fw. Generally, 
there are two possibilities for what happens when the 
fluid velocity passes the sound barrier, depending on 
whether the velocity accelerates or decelerates. If the 
fluid velocity accelerates from subsonic to supersonic, 
then information can propagate upstream from the 
sonic point, damping any tendency to develop large 
accelerations, and allowing the fluid to pass smoothly 
through the sonic point where V = iy/w. If on the other 



hand the fluid velocity decelerates from supersonic to 
subsonic, then information cannot propagate upstream, 
and in general the fluid cannot pass smoothly through 
a sonic point. Instead, the fluid steepens into a shock, 
and the velocity V changes discontinuously from being 
supersonic, \V\ > y/w, to being subsonic, \V\ < ^/w. 



G. Similarity differential equations 

Although r/t would seem to be a natural choice 
of dimensionless integration variable, in practice it 
proves to be a poor choice, because r/t does not vary 
monotonically, but rather oscillates through zero (the 
time coordinate t oscillates through ±oo), sometimes 
many times, as the baryonic fluid inside the black hole 
transitions between ingoing and outgoing. A suitable 
alternative choice of dimensionless integration variable 
is the dimensionless time parameter x defined by 



dx 



dr 
r 



(93) 



evaluated along the path of the baryonic fluid. The 
dimensionless time parameter x naturally increases 
monotonically, since the proper time r does. The 
baryonic proper time r, the time coordinate t, and the 
radial coordinate r evolve along the path of the baryonic 
fluid as 



dr 

dx 
dint 

dx 
G?lnr 

dx 



V 



(94a) 
(94b) 
(94c) 



The time coordinate t is not used in the results section, 
flIVI but the differential equation l|94b|l governing its 
evolution is given for completeness. Equation (|94b|) 
presumes that the gauge of baryonic time t is chosen 
in the natural way, such that the units of time are the 
same as the units of radius, so that r/t is a dimensionless 
variable. 

An overcomplete set of equations [only three of the 
four equations l|95() below are independent, the four 
variables £*, £ r , 7, and (3 being related by +7£ r = 1, 
equation (|97a|l below] governing the self-similar evolution 
of the remaining variables is 



dt? 
dx 
dC 
dx 
( h 
dx 
dfi_ 
dx 



72/ — (l+w)zj£ t r 



(95a) 
(95b) 
(95c) 
(95d) 
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together with 



dln[r 1+3w z(£ r ) 1+w ] 2z e s 



dx 



dlnQ 
dx 



z 
— s 



(96a) 
(96b) 



To maintain numerical precision, it is important to 
avoid expressing small quantities as differences of large 
quantities. A suitable choice of variables to integrate 
is + £ r , — 7, 7, each of which can be tiny in 
some circumstances. Starting from these variables, the 
following chain of equations yield the remaining variables 
in a fashion that ensures numerical precision: 



e - e = 



2-(/j+7)(e*+r) 



/?-7 

— = l + (/3 + 7 )(/3-7) 



(97a) 
(97b) 
(97c) 



For reference, the differential equation for X , useful in 
ray- tracing, see equation ijTUI) . is 

— = -r 



H. Boundary conditions at the outer sonic point 

The boundary conditions of the calculation are set 
in this paper at an outer boundary, taken to be a 
regular sonic point, outside the outer horizon of the 
black hole, where the infalling baryonic fluid transitions 
smoothly from subsonic to supersonic velocity. The 
behaviour in the vicinity of sonic points in general 
relativistic similarity solutions has been discussed by 

In setting the boundary conditions at the outer sonic 
point, we choose not to enquire how the accreting fluid 
managed to get to that point, but simply assume that 
astrophysical processes can arrange themselves so as to 
feed the black hole in a steady self-similar fashion. Light 
curves of real accreting black holes are typically variable 
rather than steady [72, IzHi which could indicate that 
accretion flows are typically non-steady, but we ignore 
this difficulty, assuming that at least in principle a 
black hole could accrete steadily. Integrating outwards 
from the sonic point reveals that the similarity solutions 
typically do not extend to infinity, but rather terminate, 
either at a stagnation point where the velocity V is 
zero and the fluid wants to turn back on itself, or at 
an irregular sonic point where the acceleration diverges. 
Thus the similarity solutions considered in this paper 
are typically not complete self-consistent solutions valid 
to arbitrary distances from the black hole. Again, we 
choose to ignore this difficulty, noting that the conditions 
for self-similarity, such as a relativistic equation of state, 
could well break down, so the failure of the solutions to 
extend to infinity is not necessarily a fatal difficulty. 



At the sonic point, where the fluid velocity equals 
the sound speed, V = y/w, the denominator of the 
expression (|92|l for the dimensionless acceleration y is 
zero, and the numerator must simultaneously vanish 
for the acceleration to remain finite. The vanishing of 
the numerator and denominator of the right hand side 
of equation (|92|l imposes two boundary conditions at 
the sonic point, and a third condition follows if y is 
taken to be not only continuous but also differentiable 
at the sonic point. Physically, sound waves generated 
by discontinuities near the sonic point can propagate 
upstream, modifying the flow so as to ensure a smooth 
transition through the sonic point. The value of the 
dimensionless acceleration y at the sonic point, where 
the numerator and denominator of equation l|92|) vanish, 
is given by the ratio of the derivatives of the numerator 
and denominator, according to L'Hopital's rule. If these 
derivatives are expanded according to equations (|95|) and 
(|96|l . then the result is an equation of the form 



a + by 
c + dy 



(99) 



where a, 6, c, and d are functions of dimensionless 
variables not including y. Equation l|99[) is a quadratic 
equation for y, and generically there are two solutions, 
if solutions exist. Numerically, when two solutions exist, 
one represents a transition from subsonic to supersonic, 
while the other represents a transition from supersonic 
to subsonic; it is the former that is physically relevant. 

For a sonic point at the outer boundary to be 
acceptable, it must satisfy four conditions. First, the 
sonic point must be regular, which requires that a 
solution to equation l|99|l exist. Second, the baryons must 
fall inward, that is, the velocity V of the similarity frame 
relative to the baryons must be positive, V = +-\/w. 
Third, the velocity V must transition from subsonic to 
supersonic as the baryons fall inward, which is a condition 
on the derivative of V. Fourth, the radial 4-gradient 
(/3, 7) = (<9 t r, d r r) must be spacelike, 7 2 — 1 > 0, which, 
as proven in Appendix [BJ is a necessary condition that 
the boundary be causally connected to (not separated by 
a horizon from) a hypothetical asymptotically flat empty 
region of space at large radius r. In the cases we have 
investigated, it is the fourth condition, -f 2 — 1 > 0, that 
sets limits on the values of physical parameters at the 
outer sonic point. 

Given the assumption that the outer boundary is 
a regular sonic point, there are two further boundary 
conditions to be set for the baryons: the dimensionless 
accretion rate rj s = ar/t, and the charge-to-mass ratio 
Q/M c of the black hole, where M c is the charge- 
augmented mass defined by equation (|100|l below. 

The dimensionless accretion rate rj s = ar/t at the sonic 
point is roughly speaking the velocity at which the black 
hole is expanding, the coefficient a being a constant of 
order unity whose value depends on how one chooses to 
fix the gauge of time t. For most calculations in this 
paper and in Paper 2 we choose a rather large value of the 
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FIG. 1: Charge-to-mass Q/M c of the black hole as a function 
of the charge-to-mass density q/pb of baryons at the outer 
sonic point boundary, for an accretion rate rj s = 0.1. Two 
cases are shown, one with zero conductivity, k = 0, the other 
with finite conductivity, k = 1. The charge-to- mass Q/M c 
of the black hole generally increases as the charge-to-mass 
density q/pt at the outer boundary is increased, although 
in the zero conductivity case the charge-to-mass Q/M c does 
turn down slightly at the largest charge-to-mass densities 
q/pb- The maximum charge-to-mass density q/pb is set by 
the condition that 7 2 — 1 > 0. 

accretion rate, r/ s = 0.1, which corresponds roughly to a 
black hole expanding at a tenth of the speed of light . In 
astronomically realistic black holes, where the accretion 
timescale is long compared to the characteristic timescale 
of the black hole, the accretion rate T] s is unlikely to be 
this large. In the calculations reported in Ej iUVHVI we 
choose to use an unrealistically large value of r/ s partly 
so as to amplify the difference between the similarity 
solutions and the vacuum black hole, and partly to avoid 
the risk of numerical artifacts that might possibly be 
associated with a tiny r) s . 

The charge-augmented mass M c at radius r is defined 

by 

M c = M + ^— (100) 

where M and Q are the mass and charge interior to radius 
r. In the Reissner-Nordstrom geometry, the charge- 
augmented mass M c is constant as a function of radius. 
The mass Q 2 /2r = J r ° ° \(Q /r 2 ) 2 /8w] Aitr 2 dr on the right 
hand side of equation ((100(1 is the mass-energy that would 
be in the electric field Q /r 2 outside radius r if there were 
no charge outside r. The charge-augmented mass M c at 
the sonic point is the mass that the black hole would 
appear to have at infinite distance if there were no mass 
or charge outside r. In this paper, values of radius r are 
reported in units where M c — 1 at the outer sonic point. 
This is not a boundary condition on M c , just a choice of 
units. 

The charge-to-mass Q/M c of the black hole at the 



outer sonic point depends on the entire accretion history 
of the black hole. In the similarity solutions, this 
charge-to-mass Q/M c is determined by the charge-to- 
mass density q/ Pb of the baryons being accreted into 
the black hole. Figure ^ shows that the charge-to-mass 
Q/M c of the black hole at the outer sonic point generally 
increases as the charge-to- mass density q/ pb of baryons 
at the outer sonic point is increased, although Q/M c does 
turn down slightly at the largest values of q/pb when the 
conductivity is small. As the conductivity of the baryonic 
plasma is increased, a larger charge-to-mass density q/pb 
at the sonic point is needed to produce a given charge- 
to-mass Q/M c of the black hole. This is because a larger 
conductivity allows the charge-to- mass density q/pb that 
is actually accreted through the outer horizon of the black 
hole to be significantly smaller than the charge-to-mass 
density at the sonic point. 



IV. RESULTS 

This section presents results for black holes that 
accrete only baryons. In Paper 2 [37J the black hole is 
allowed to accrete dark matter in addition to baryons. 

Geometric units G = c — M c = 1 are used here 
and throughout the remainder of this paper, where M c , 
equation 1 (100(1 , is the charge-augmented interior mass of 
the black hole evaluated at the outer boundary, the sonic 
point. 



A. Uncharged black hole 

Figure |21 shows the simplest of the black hole solutions 
considered in this paper, that for a black hole which 
accretes baryons with zero charge. Accreting no charge, 
the black hole is, naturally, uncharged. 

The uncharged solutions are characterized by a single 
free parameter, the dimensionless accretion rate rj s at the 
outer boundary, the outer sonic point. As with all the 
models illustrated in this paper, the accretion rate rj s is 
set equal to 0.1 at the outer sonic point 

Vs = 0.1 (101) 

which roughly speaking means that the black hole is 
expanding at a tenth of the speed of light. This accretion 
rate is large compared to that of an astronomically 
realistic black hole, but the large value makes it easier to 
discern the difference in geometry between the similarity 
solutions and a vacuum black hole. Results for smaller 
accretion rates are qualitatively similar. The accretion 
rate of the uncharged black hole is limited to i] s < 0.3616 
by the constraint j 2 — (3 2 > at the outer sonic point. 

Tabic [I] compares the radii r of several points in 
the similarity solution to those in the Schwarzschild 
geometry. As can be seen, there is a considerable degree 
of commonality, despite the 'large' accretion rate rj s = 0.1 
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FIG. 2: (Color online) The simplest model considered in 
this paper: a black hole that accretes uncharged baryons, 
at rate r\ a = 0.1. Quantities are plotted against radius r in 
units where the charge-augmented mass at the outer sonic 
point is unity, M c = 1. Lines are dashed where quantities 
are negative. Disks mark the outer sonic point, where V = 
\/w « 0.577, at which the boundary conditions are set. Short 
horizontal bars mark the horizon, where V = 1. (Upper 
panel) Proper velocity V of the similarity frame relative to 
the baryonic frame. (Middle panel) Dimensionless proper 
baryonic mass density z = AiYr 2 pb. (Bottom panel) Interior 
mass M, proper acceleration g experienced by the baryonic 
fluid, and the homothetic scalar H, equation l|58|l . 



in the similarity solution. In the similarity solution, the 
radii are as measured in a frame falling with the baryons 
(the radius r is not a similarity variable, hence the need 
to specify the reference frame). In the Schwarzschild 
geometry, the position of the outer sonic point is given 
for a tracer relativistic fluid (p/p = 1/3) of uncharged 
particles which free-fall from zero velocity at infinity. 



TABLE I: Radii of events in an uncharged black hole 
accreting baryons at rate rj s =0.1, compared to corresponding 
radii in the Schwarzschild (Schw) solution. 



r = 




FIG. 3: (Color online) Penrose diagram for the uncharged 
black hole of Figure The arrowed line represents 
schematically the worldline of a parcel of baryonic fluid. The 
thin grey lines are lines of constant self-similar coordinate. 
The self-similar solution extends some way outside the sonic 
point outside the outer horizon, but it does not continue to 
r — oo; the Penrose diagram as drawn presumes that the 
spacetime joins to an asymptotically flat space as r — > oo. 



As illustrated in Figure [21 with no charge to repel 
its fall, the baryonic fluid inside the black hole plunges 
straight to a singularity at zero radius, the infall velocity 
diverging as V ~ r _1 as r — > 0. As the fluid falls, the 
dimensionless proper baryonic mass density z = Aitr 2 pt> 
tends to a constant, z — > 0.02092, indicating that the 
baryonic density diverges as pi, ~ r~ 2 as r — > 0. 

The baryons develop a pressure gradient which causes 
a mild outward proper acceleration g, but this is not 
enough to prevent gravitational collapse. It takes a 
proper time of At = 2.164 for the baryons to fall from 
the horizon to the singularity, somewhat larger than the 
proper time At = 4/3x 1.333 for a test particle that free 
falls radially from zero velocity at infinity to fall from the 
horizon to the singularity of a Schwarzschild black hole. 

Among the variables shown in Figure |21 is the 
homothetic scalar H, equation (|58|l . Besides having 
the virtue of being, like the interior mass M, a 
gauge-invariant scalar, the homothetic scalar H plays 
a fundamental role in describing geodesies, see flIII Dl 
and hence in determining what things look like if you 
fall into a black hole, a question addressed in f|V] In 
particular, the place where the homothetic scalar H 
reaches a maximum sets the location of the equivalent 
of the photon sphere, which determines the angular size 
of the black hole perceived by an infalling observer. 
Although the boundary conditions are set at the outer 
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sonic point, Figure [5] also shows a partial continuation 
of the similarity solution outward from the outer sonic 
point. 

Figure [3| shows a Penrose diagram of the uncharged 
black hole. The Penrose diagram looks qualitatively 
similar to that of the Schwarzschild solution, except 
that whereas in the Schwarzschild solution the lower 
left diagonal edge of the Penrose diagram marks the 
antihorizon, in the similarity solution the same diagonal 
edge marks the collapse event at zero radius. The 
similarity solution does not really include the collapse 
event itself, and presumably the Penrose diagram would 
be changed if the spacetime incorporated the collapse 
event. 



B. Black hole accreting charged, non-conducting 
baryons 

Figure 0] shows a similarity solution for a black hole 
which accretes charged baryons with zero conductivity. 
There are two free parameters, the accretion rate 
77s, which is set to the same value 0.1 as before, 
equation (|101fl . and the charge-to- mass ratio Q/M c of 
the black hole, which is set at the outer sonic point to 



Mr 



0.8 



(102) 



where M c is the charge-augmented mass defined by 
equation 1)100(1 . The condition j 2 — 1 > at the outer 
sonic point limits the maximum possible accretion rate 
and maximum charge-to- mass density at the sonic point. 
At fixed Q/M c = 0.8, the accretion rate is limited to 
T] s < 0.3496. At fixed r/ s = 0.1, the charge-to-mass 
ratio is limited to Q/M c < 1.0060, the maximum value 
being attained at slightly less than the maximum charge- 
to-mass density q/pb, whereat Q/M c 0.99979536, as 
illustrated in Figure ^ 

The black hole charge Q is taken without loss of 
generality to be positive. The charge of real black holes is 
most likely positive, albeit minuscule, because the larger 
mass-to-charge of protons compared to electrons makes 
it easier for a black hole to accrete positive charge. 

Unlike the uncharged black hole, the baryons inside the 
charged black hole do not fall to a singularity. Instead, 
they are repelled by the charge of the black hole, become 
outgoing, and drop through the Cauchy horizon. The 
transition from ingoing to outgoing is marked by the 
proper velocity V of the similarity frame relative to the 
baryons passing from +oo to — oo. As described at the 
end of fflll CI the fluid remains perfectly well-behaved 
through this point: the velocity V passes through infinity 
because the homothetic 4- vector £ in the baryonic frame 
switches from pointing forwards to backwards in time. 

Table [n] compares the radii r of several points in the 
similarity solution to those in the Reissner-Nordstrom 
(RN) geometry with the same charge-to-mass 0.8. As 
previously found in the uncharged model of illV Al there 
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FIG. 4: (Color online) Similar to Figure [5] but for a black 
hole with charge-to- mass Q/M c — 0.8, accreting baryons with 
zero conductivity, k = 0. Lines are dashed where quantities 
are negative. Disks mark sonic points, where V = ±\/w. 
Short horizontal bars mark horizons, where V = ±1. (Upper 
panel) Proper velocity V of the similarity frame relative to 
the baryonic frame. The velocity passes from V = +oo 
to — oo as the charged baryonic fluid passes from ingoing 
to outgoing. The baryonic fluid drops through the Cauchy 
horizon, beyond which the solution terminates at an irregular 
sonic point. (Middle panel) Dimensionless proper baryonic 
mass and charge densities z = 4Tvr 2 pb and z q = 4-nr 2 q. 
(Bottom panel) Interior mass M, proper acceleration g, and 
the homothetic scalar H. 



is a considerable degree of commonality, despite the 
'large' accretion rate ?7 S = 0.1 in the similarity solution. 
In the similarity model, the radii are as measured in a 
frame falling with the baryons. In the RN geometry, the 
positions of sonic points, and the places where charged 
baryons become outgoing, and where 7 goes negative, 
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are given for a tracer relativistic fluid (p/p = 1/3) of 
particles whose charge-to-mass is the same as that of the 
black hole, 0.8, and which free-fall from zero velocity at 
infinity. 

TABLE II: Radii of events for the model with accretion rate 
r/s — 0.1, charge-to-mass Q/M c = 0.8, and zero conductivity, 
compared to corresponding radii in the Reissner-Nordstrom 
solution with the same charge-to-mass Q/M c — 0.8. 



Event 


Model 


RN 


Photon sphere 


2.439 


2.485 


Sonic point 


2.439 


2.122 


Outer horizon 


1.813 


1.6 


(/3, 7) becomes timelike 


1.643 


1.6 


Charged particles become outgoing 


0.639 


0.64 


7 becomes negative 


0.614 


0.64 


Inner horizon 


0.3782 


0.4 


(/3,7) becomes spacelike 




0.4 


Sonic point 


0.3775 


0.3994 



If the baryons could see ingoing light or matter falling 
from the outside universe, then as the baryons passed 
through the Cauchy horizon they would see the light 
or matter infinitely blucshifted. However, it is being 
assumed here that there is no ingoing radiation or 
matter. Physically, one can imagine that the baryons are 
sufficiently opaque that no ingoing radiation or matter 
reaches the baryons near the Cauchy horizon. In Paper 2 
the effects of ingoing dark matter, whose streaming 
through the outgoing baryons leads to mass inflation near 
the Cauchy horizon, will be considered explicitly. 

Inside the Cauchy horizon the baryonic fluid continues 
to decelerate until the flow velocity V equals the speed 
of sound, at which point the proper acceleration g of the 
fluid diverges. Normally this is a signal that a shock must 
form, which decelerates the fluid discontinuously from 
supersonic to subsonic. The position of the shock is a free 
parameter, which in practice is constrained over a rather 
narrow range of radii where the fluid is subluminal but 
supersonic. Normally the position of the shock would be 
fixed by requiring that the solution downstream continue 
in a regular fashion. In the case under consideration, 
the gas, having been decelerated below the speed of 
sound by the shock, subsequently accelerates back up 
to the speed of sound, and normally one would require 
that the place where the velocity accelerates back from 
subsonic to supersonic be a regular sonic point, where the 
acceleration is finite and preferably differentiable. Here 
however the acceleration diverges for all values of the 
shock position, so there is no consistent continuation of 
the similarity solution. 

Figure \E\ illustrates three possible shocks, ranging 
from extremely weak (small shock velocity) to extremely 
strong (highly relativistic shock velocity). In all cases, 
whereas the pre-shock fluid decelerates strongly, the 
post-shock fluid re-accelerates strongly inward, and soon 
accelerates back up to an irregular sonic point, where the 
acceleration diverges. 
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FIG. 5: (Color online) Putative continuations of the 
similarity solution inside the Cauchy horizon for the model 
shown in Figure [I] which has Q/M c — 0.8 and zero 
conductivity. Three cases are shown, ranging the gamut 
from an extremely weak shock to an extremely strong shock. 
In all cases the shocked fluid continues downstream for a 
short distance before terminating at an irregular sonic point, 
marked by a disk, where the proper acceleration diverges and 
the similarity solution cannot be continued. Dotted lines 
connect pre-shock and post-shock quantities. (Top panel) 
Proper velocity V of the similarity frame relative to the 
baryonic frame. The Cauchy horizon, where V = —1, is 
marked by a short horizontal bar. A shock decelerates the 
baryonic fluid from supersonic to subsonic, through the sound 
speed where V — — \/w ~ 0.577. (Bottom panel) Proper 
acceleration g experienced by the baryonic fluid in the vicinity 
of each of the three putative shocks. The acceleration changes 
from outward (negative, dashed lines) in the pre-shock fluid 
to inward (positive, continuous lines) in the post-shock fluid. 



The failure of the similarity solution to continue to zero 
radius is discussed in i llVDI 



C. The radial 4-gradient is timelike inside the 
Cauchy horizon 

An apparently small but nevertheless crucial difference 
in Table |H] between the similarity solution and the 
RN geometry concerns the radial 4-gradient (/3, 7) = 
(dtr,d r r). In the RN geometry the radial 4-gradient, 
having switched from spacelike to timelike at the outer 
horizon, then switches back to spacelike at the inner 
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horizon. By contrast, in the similarity solution, the radial 
4-gradient, having switched from spacelike to timelike a 
little way inside the outer horizon, never changes back to 
being spacelike. 

The timelike character of the radial 4-gradient ((3, 7) 
in the similarity solution means that it is impossible 
inside the Cauchy horizon to accelerate (with rockets, 
say) to a frame which is at rest in radius r, with (3 — 0: 
all locally inertial frames inside the Cauchy horizon are 
falling radially inward, (3 < 0. Thus a person inside the 
Cauchy horizon is doomed to move inward to smaller 
radius r, which is quite unlike the RN solution. Because 
the similarity solutions do not continue to zero radius, it 
is not possible to say exactly what happens, but if {{3, 7) 
remains timelike down to zero radius, then an infaller will 
hit zero radius, r = 0, in a finite proper time. The best 
that an infaller can do to delay the inevitable fate is to 
accelerate to a frame where 7 = 0, the frame where (3 is 
least negative. 

This seems paradoxical: if the region inside the Cauchy 
horizon is subluminal, why can't an infaller stay away 
from zero radius? In the context of similarity solutions, 
subluminal means being able to move either outward or 
inward relative to the similarity frame. Here however the 
similarity frame is itself contracting inside the Cauchy 
horizon. The similarity frame moves with r oc t, 
and normally one would think that this means that 
the similarity frame is expanding. However, inside the 
Cauchy horizon it is consistent to think of the coordinate 
time t as being negative, and increasing towards t = 0; 
alternatively, one can think of the coordinate time t 
as being positive and decreasing towards t — (the 
ambiguity in the sign of t expresses gauge freedom; 
the gauge-invariant statement is that 77 = ar/t must 
be negative). The radius r of the similarity frame 
thus contracts to r = at t — 0. In the ligh t of 
the Penrose-Hawking singularity theorems [35t H3L l65| 
one might reasonably expect a presumably spacelike 
singularity at zero radius, although this expectation 
cannot be confirmed here explicitly, because of the failure 
of the similarity solutions to continue consistently to zero 
radius. 

Parcels of baryonic fluid that fall through the Cauchy 
horizon later find themselves at larger radius r. Thus, 
again paradoxically, even though the similarity frame is 
contracting from the point of view of objects at rest 
in the similarity frame inside the Cauchy horizon, the 
Cauchy horizon is nonetheless expanding in the sense 
that baryons that fall in later find the Cauchy horizon 
at larger radius r. 

Figure [B] shows a tentative Penrose diagram of the 
spacetime under consideration. It is tentative because 
it extrapolates to a spacelike central singularity, which is 
not established by the similarity solution. The Penrose 
diagram presumes that the baryons undergo a self- 
similar shock before falling to the central singularity. 
A question mark on the diagram emphasizes the fact 
that the similarity solution leaves undetermined how the 
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FIG. 6: (Color online) Tentative Penrose diagram for a 
charged black hole, such as that shown in Figures and 
|5J in which the baryons drop through a Cauchy horizon 
and undergo a shock before falling to a central spacelike 
singularity. The arrowed line represents schematically the 
worldline of a parcel of baryonic fluid. The thin grey lines 
are lines of constant self-similar coordinate. The self-similar 
solution extends some way outside the sonic point outside 
the outer horizon, but it does not continue to r = 00. The 
Penrose diagram as drawn presumes that the spacetime joins 
to an asymptotically flat space as r — > 00. Similarly, the self- 
similar solution extends some way inside the shock inside the 
Cauchy horizon, but it does not continue to a singularity, but 
rather terminates at an irregular sonic point. The Penrose 
diagram as drawn presumes that the spacetime terminates at 
a spacelike singularity at r = 0. A question mark emphasizes 
the fact that the similarity solution leaves undetermined how 
the spacetime continues (or not) to the singularity. 



spacetime continues (or not) to the singularity. The 
diagram shows that if a self-similar shock is present inside 
the Cauchy horizon, then the shock propagates inwards 
from parcels of baryonic fluid accreted at later times 
towards parcels of baryonic fluid accreted at earlier times. 

The left (lower and upper) diagonal edges of the 
Penrose diagram mark the collapse event at r — 0. The 
similarity solution does not really include the collapse 
event itself, and presumably the left part of the Penrose 
diagram would be changed if the spacetime incorporated 
a realistic collapse event. Notice that, according to the 
Penrose diagram, a test particle that falls into the black 
hole and remains ingoing will encounter a null singularity 
at zero radius. It is a basic hypothesis of the similarity 
solution that such ingoing particles are not actually 
present: if any finite energy density of ingoing particles 
were present, then that would lead to mass inflation near 
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the inner horizon, as considered in Paper 2. Nevertheless, 
if the similarity solution is supposed to remain valid all 
the way to the collapse event at r = 0, then an ingoing 
test particle with a zero energy-momentum tensor would 
in principle approach r = 0, and in so doing would pass 
through outgoing baryonic fluid accreted at ever closer 
to the initial collapse event. 

Question for the reader: If you fall inside the Cauchy 
horizon, you can hover just below the Cauchy horizon 
by accelerating hard — but in which direction? Before 
answering this question let us orient ourselves inside the 
black hole. Before you fall into the black hole, equip 
yourself with a gyroscope that is initially set to point 
radially inward into the black hole. As you fall inward, 
you define the direction in which the gyroscope points 
to be the immutable direction towards the black hole. 
The answer to the question posed is then: In order to 
hover just below the Cauchy horizon, you must accelerate 
inwards towards the black hole, in the direction indicated 
by the gyroscope. You might have thought that you 
would have to accelerate outwards, in the direction of 
baryons that fall through the Cauchy horizon after you, 
but this is incorrect: you must accelerate inwards, in the 
direction of baryons that fell in before you. Even though 
from your own point of view you are hovering at the 
Cauchy horizon, a person who falls through the Cauchy 
horizon after you never sees you, either at the Cauchy 
horizon or anywhere else. On the other hand a person 
who falls through the Cauchy horizon before you does 
see you, not necessarily close to the Cauchy horizon from 
their point of view, rushing rashly by in the direction 
towards the black hole. You can draw these conclusions 
by examining the Penrose diagram in Figure El 

It is worth emphasizing how subtle is the small 
difference between the similarity and RN solutions that 
leads to so dramatically different causal behavior inside 
the Cauchy horizon. Whereas in the RN solution the 
radial 4-gradient (/3, 7) changes from being timelike to 
spacelike at the Cauchy horizon, in the similarity solution 
(/?, 7) remains timelike at and inside the Cauchy horizon. 
Since 2M/r— 1 = (3 2 — 7 2 , equation J2U, the timelike or 
spacelike character of (/3, 7) depends on whether M/r is 
greater or less than 1/2. Whereas in the RN solution 
M/r = 1/2 at the Cauchy horizon, in the similarity 
solution M/r = 0.510 at the Cauchy horizon, just 
slightly greater than 1/2. If the similarity solution is 
continued inside the Cauchy horizon, then M/r > 0.508 
at the terminal sonic point, the lower limit M/r = 
0.508 occurring for the case where there is no shock. 
Unfortunately, a proof that (/3, 7) must remain timclike 
inside the Cauchy horizon, along the lines of that in 
Appendix [51 fails; but numerically ((3, 7) does remain 
timelike, and if it does so all the way to zero radius, 
then objects that fall through the Cauchy horizon must 
inevitably fall to zero radius. 

We have experimented with models in which the black 
hole charge-to- mass Q/M c is at and near its maximum 
value, but the models do not differ qualitatively from 



that shown in Figure 01 The charge-to-mass is limited 
by the constraint j 2 — (3 2 > 0, which is the condition, 
Appendix El that the outer sonic point not be separated 
by a horizon from a hypothetical asymptotically flat 
empty region of space at large radius r. The constraint 
intervenes before the black hole becomes extremal, that 
is, before the inner and outer horizons of the black 
hole become coincident (which happens in the Reissner- 
Nordstrom geometry when Q/M c = 1). There is 
no sign of phenomena associated with cricital collapse 
j2^1 123 . Isfj . such as ringing, or the appearance of a naked 
singularity. 



D. Incompleteness of solutions inside the Cauchy 
horizons 

In 31 V Bl it was found that the similarity solution inside 
the Cauchy horizon could not be continued consistently 
to zero radius, but rather terminated at an irregular sonic 
point, where the acceleration diverged. The solution 
terminated whether or not a shock was introduced. We 
find the same behavior whenever the similarity solution 
drops inside the Cauchy horizon: in all the cases that 
we have examined, including those described in the rest 
of this paper, and in §IVB of Paper 2, the solution 
terminates at an irregular sonic point, with or without a 
shock, just inside the Cauchy horizon. 

Does this mean that the entire similarity solution must 
be discarded as inconsistent? The methods considered in 
this paper are insufficient to supply a definitive answer 
to this question, but physically it seems reasonable 
that the similarity solution outside the Cauchy horizon 
should be valid even if there is no consistent self-similar 
continuation inside the horizon. Since no information 
can propagate outward from the Cauchy horizon, the 
solution outside the Cauchy horizon cannot know that 
the solution fails inside the Cauchy horizon. 

As remarked above, the similarity solution does not 
include the instant of gravitational collapse where the 
black hole first forms, at M = r = t = 0. But if a seed 
black hole, having formed with some geometry or other, 
accretes in a self-similar fashion, then it seems reasonable 
that the black hole could, after many doublings of its 
mass, settle asymptotically to the self-similar form. The 
interior mass M and charge Q of the black hole are 
generated self-consistently by the accretion of charged 
baryons, so it seems reasonable to expect that after a 
sufficiently long time the black hole would forget what 
happened at its formation. 

The presence of any outgoing tail of radiation produced 
by the gravitational collapse of the black hole pa. I2H ItH 
68) is not relevant to the question of the consistency of 
the similarity solutions being considered here. In the 
similarity solutions considered here, the baryonic fluid is 
already outgoing, and the presence of additional outgoing 
radiation cannot prevent the baryonic fluid from passing 
through the outgoing Cauchy horizon. What can prevent 



20 




FIG. 7: Velocity V versus radius r for pure bayonic models with various conductivities k, increasing from zero at top left 
to a large value somewhat less than the maximum conductivity K max ~ 67.0 at bottom right. At conductivities less than the 
continuum threshold Koo, the top six panels, the baryonic fluid drops through the Cauchy horizon unless the conductivity equals 
one of an infinite spectrum of discrete values K n , in which case the fluid falls to a singularity at zero radius. The first two cases 
of the discrete spectrum, ki and K2, are shown. The velocity V passes through ±oo where the baryonic fluid transitions between 
ingoing (V > 0, continuous lines) and outgoing (V < 0, dashed lines). Solutions between values of the discrete spectrum k„ are 
characterized by the number of times that the baryonic fluid cycles between ingoing and outgoing. At the continuum threshold 
Koo, the bottom left panel, the solution exhibits discrete self-similarity, the baryonic fluid cycling repeatedly between ingoing 
and outgoing. At and above the continuum threshold k^, the bottom three panels, all solutions fall to a central singularity. 
Disks mark sonic points, where V = ±\/w m ±0.577. Horizons occur where V = ±1. 



the baryonic fluid from dropping through the Cauchy 
horizon is ingoing matter or radiation (which leads to 
mass inflation; see Paper 2). To have a causal effect on 
a parcel of outgoing baryonic fluid, such ingoing matter 
must necessarily be accreted after the baryons. It is thus 
evident that what happened to the black hole at its birth 
cannot affect whether or not the baryonic fluid drops 
through a Cauchy horizon, though the birth event clearly 
can affect what happens in the subluminal region inside 



the Cauchy horizon. 



Whatever the case, the results of the present paper 
provide a definite set of self-similar boundary conditions 
at the Cauchy horizon, which could be supplied to a 
general relativistic hydrodynamic code. Here we leave 
the intriguing question of what happens beyond the 
Cauchy horizon to a future investigation. 
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the value given for K max is theoretical (inferred from the 
condition 7 2 — /3 2 = 0); numerical integration becomes 
unstable if n is within about 0.1 of K max sa 67.0. 



TABLE III: Spectrum of critical conductivities for a black 
hole with accretion rate rj s =0.1 and charge-to-mass Q/M c = 
0.8. 



K.2 
K3 



0.9982939540115 
1.1177942541104 
1.1654295843302 
1.1903424212409 
1.30382 

66.998112873734 



FIG. 8: (Color online) Velocity times radius, Vr, for two 
models with conductivities k = ki ± e just either side of 
the value rei ~ 0.998 that separates solutions that cycle once 
versus twice between ingoing (V > 0, continuous lines) and 
outgoing (V < 0, dashed lines). The inset shows a detail of 
how the velocities in the two models peel off from the critical 
solution, which itself would continue on the outgoing path as 
a horizontal line to zero radius. 



E. Black hole accreting charged, conducting 
baryons 

Figure [7] shows the proper velocity V of the similarity 
frame relative to the baryonic frame for black holes ac- 
creting baryons with a range of electrical conductivities. 
The accretion rate t] s — 0.1 and black hole charge-to- 
mass Q/M c = 0.8 are set to the same values as before, 
equations (|101|l and (|102|l . 

The solutions reveal a surprisingly rich structure as the 
dimensionless conductivity coefficient k, equation l|48|l . is 
increased. The baryonic fluid accreted by the black hole 
undergoes one of two possible fates: either the baryonic 
fluid drops through a Cauchy horizon, or else it falls to 
a singularity at zero radius. Mostly, the former fate 
happens if the conductivity is less than a continuum 
threshold conductivity Koo « 1.304, while the latter fate 
happens if the conductivity is greater than the threshold. 
However, the baryons fall to a central singularity not 
only in the continuous range of conductivities > Koo , but 
also at any of an infinite spectrum of discrete values K n 
(n = 1,2,...), asymptoting to the continuum threshold 

Table ITTTl lists the first four of the discrete values k„, 
the continuum threshold value Kqo, and the maximum 
conductivity K max . The maximum conductivity K max 
is set by the requirement that 7 2 — /3 2 > at the 
sonic point, which is the condition that the sonic point 
be causally connected to (not separated by a horizon 
from) a hypothetical asymptotically flat empty region 
of space at large radius (see Appendix [BJ). Whereas 
the values for K n and in Table IIIII are determined 
numerically (with fewer significant digits for because 
of the greater difficulty in locating this point precisely), 



Solutions between the discrete values Kyi &X6 charac- 
terized by the number of times n that the baryonic 
fluid transitions between ingoing and outgoing. At 
conductivities below the first critical value, k < K\, 
the fluid transitions from ingoing to outgoing and then 
drops through the Cauchy horizon. At conductivities 
between the first two critical values, k\ < k < K2, the 
fluid transitions from ingoing to outgoing to ingoing to 
outgoing before dropping through the Cauchy horizon. 

Collapse to a singularity at one of the discrete 
conductivities n n requires extreme fine-tuning: the 
conductivity must be exactly equal to the critical value. 
With any tiny departure from the critical conductivity, 
the baryons will not fall to a singularity, but will instead 
drop through the Cauchy horizon. 

Figure [S] shows two solutions whose conductivities are 
just slightly above and below, by a few x 10~ 14 , the first 
critical value K\. The two solutions initially track each 
other closely. The velocity V passes through ±oo as 
the baryons switch from ingoing to outgoing, and then 
settles to a power law V cx r . If the conductivity 
were exactly critical, k = k%, then the velocity would 
continue as V oc down to zero radius. In the 

solutions shown, the velocities depart from the power 
law at a small radius r ~ 10~ 15 . There is nothing 
special about this radius: the departure occurs at smaller 
and smaller radius as k approaches K\. In the solution 
just below k±, once the velocity V departs from the 
power law, it decreases rapidly, and the baryons soon 
drop through the Cauchy horizon. Conversely, in the 
solution just above K\, once the velocity V departs from 
the power law, it increases rapidly, passes from — oo to 
+oo as the baryons transition back from outgoing to 
ingoing, and settles back to another power law V oc r . 
The velocity remains on this power law for some time 
before steepening into another power law V oc r~ 4 , then 
steepening again to pass from +oo to — oo as the baryons 
transition for a second time from ingoing to outgoing. 
The velocity then decreases rapidly, and the baryons 
finally drop through the Cauchy horizon, at a radius some 
35.4 decades smaller than the radius at which the below- 
critical solution already dropped through. 

The behavior of the solutions close to the critical 
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FIG. 9: (Color online) Similar to FiguresHandQ] Model with 
conductivity k = 0.8, less than the first critical conductivity 
«i. The baryonic fluid drops through the Cauchy horizon. 



conductivity n± is similar to that nicely described by and 
illustrated in Figure 1 of Gundlach (2003) [3(| in the 
context of critical gravitational collapse. 

Frolov and Pen (2003) j^D] remark that the solutions 
just above and below the critical solution in critical 
gravitational collapse appear quite different. They pose 
this as a paradox, and conjecture that perhaps the two 
solutions differ only by a gauge change. The results 
obtained here suggest a different interpretation, that two 
spacetimes either side of a critical point may indeed be 
radically different. As the critical point is approached, 
the differences, though still present, recede to ever 
smaller radius. 

It seems possible that the oscillatory character of 
the solutions with conductivities between k\ and 
might perhaps be related in some way to the oscillatory 



character of the curvature near the inner horizon of a 
rotating black hole reported by [63| . 

Figure illustrates a model with conductivity k = 0.8 
below the first critical threshold k\. Like the zero 
conductivity model shown in Figure 21 the baryonic fluid 
ends up dropping through the Cauchy horizon. The 
finite conductivity allows currents to flow so that, as 
Figure ED shows, the interior charge Q decreases along 
the path of the infalling baryons, whereas in the non- 
conducting case the interior charge Q would be constant. 
The conductivity is however not so large as to allow the 
baryonic fluid to neutralize itself entirely, so Q remains 
positive, and, like the non-conducting case shown in 
Figure the baryons end up dropping through the 
Cauchy horizon. 

In the conducting model shown in Figure the 
dimensionless charge density z q is negative inside a 
certain radius, whereas in the non-conducting model 
of Figure the charge density remained everywhere 
positive. How can the interior charge Q be positive if 
the charge density inside a certain radius is everywhere 
negative? Because conservation of charge in the 
similarity frame depends on the current j as well as the 
charge density q. 



F. Discrete self-similarity at threshold 

The self-similar solution at the continuum threshold 
conductivity Koc, is of special interest. The solution 
displays a discrete self-similarity strikingly reminiscent 
of the discrete self-similarity in critical gravitational 
collapse first discovered numerically by Choptuik (1993) 
[23]] ; see [3^ for a review. 

Figure ^| shows several variables in the similarity 
solution at the continuum threshold /-too (or at least as 
close to the continuum threshold as we could discern, 
numerically). All the variables oscillate in phase, the 
solution at each period appearing to be a scaled version 
of the solution at any other period. 

Figure 1771 shows the velocity V scaled by r 1 - 7081 . The 
scaled velocity plotted in Figure ^] passes through 3f 
cycles, each 3.9463 decades wide in radius, covering in 
all about f 20 decades of radius. 

The discrete self-similarity is not exact. For example, 
the first cycle, which starts outside the horizon, starts 
out noticeably different in shape from the later cycles. 
Presumably the discrete self-similarity operates asymp- 
totically, the shape of the cycle converging asymptotically 
to a certain form as the number of cycles increases. 

Koike, Hara & Adachi (1995) [5l| explained the dis- 
crete self-similarity observed in gravitational collapse as 
resulting from a dominant unstable mode of perturbation 
to an unperturbed self-similar solution. This idea has 
been explored and built upon by several authors 0, 0, 
EM I52I l55| ; see f° r a review. Possibly a similar 
explanation could apply in the present case. 
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FIG. 10: (Color online) Similar to FiguresElH Model with 
conductivity at the continuum threshold conductivity k = 
Koo- The model displays discrete self-similarity. 



V. APPEARANCE OF THE BLACK HOLE 

What does it actually look like if you fall inside one 
of the black holes described in this paper? A full answer 
to this fascinating question is beyond our scope, but one 
can begin to get an idea by answering a simpler question: 
how big does the black hole appear on the sky as you 
fall into it? Or at least, how big would the black hole 
appear to you if you could actually see it through the fog 
of accreting baryons, an ominous black disk silhouetted 
against the starry universe outside? 

It would of course have to be you, not a remote 
instrument, that observed the inside of the black hole, 
since an instrument inside the horizon would not be able 
to relay its observations to you outside the horizon. You 
should be sure to fall into a supermassive black hole, 




1(T 6 1(H IO" 2 10° 10 2 
Periodically folded radius r 



FIG. 11: (Color online) Discrete self-similarity at the 
continuum threshold conductivity k, = (too. The horizontal 
axis is the radius r folded over two periods, a period being 
3.9463 decades. The vertical axis is the velocity V multiplied 
by a power of radius, r B1 . The solution as shown passes 
through 31 periods, covering a total of about 120 decades of 
radius. 



since you could easily survive the tidal forces well inside 
the black hole, whereas a stellar-massed black hole would 
dismember you long before you reached the horizon. A 
useful rule to remember is that, for an observer in circular 
orbit, the tidal force is approximately one gee per meter 
at a radius where the orbital period is one second, a fact 
which the interested reader may care to confirm. 

In a Reissner- Nordstrom black hole the perceived edge 
of the black hole is set by photons orbiting at the photon 
sphere. Such photons fall from the outside universe into 
the unstable orbit at the photon sphere, circulate the 
black hole an arbitrarily large number of times, and then 
either escape back out to infinity or else fall into the black 
hole. 

In the self-similar black holes of the present paper, 
there is no photon sphere as such, no place where photons 
can orbit in circles for ever. There is however a separatrix 
between photons from the outside universe that do or 
do not fall into the black hole, and it is this separatrix, 
which we call the photon sphere equivalent, that defines 
the perceived edge of the black hole as viewed on the sky. 

Figure ^] shows, for the three models illustrated in 
Figures [21 01 an d E3 the angular size Xph and blueshift 
of photons from the edge of the black hole, as observed 
either in the baryonic rest frame, or in a rest frame 
that free-falls radially from zero velocity at infinity (in 
the latter case, the frame is the same as that of the 
pressureless dark matter considered in Paper 2). The 
two points of view are related by a radial Lorentz boost. 
The observed angular size Xph of the black hole (the 
subscript ph signifying photons from the photon sphere 
equivalent) is given by equation (|73|l . and the observed 
blueshift of photons at the edge of the black hole is given 
by equation (|74"|> . 

The appearance of an uncharged, non-conducting 
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FIG. 12: (Color online) Angular size Xph of the black hole, equation 1731 . and the blueshift of photons at the edge of the black 
hole, equation 17411 . perceived by an observer in either the baryonic frame or a radial free-fall frame, for three models. The 
models shown are: (left) the uncharged baryonic model from Figure [5] (middle) the charged non-conducting baryonic model 
from Figure [I] (right) the discretely self-similar charged baryonic model with k = k^, from Figure [TUl Disks mark the outer 
sonic point, where V = \/w ~ 0.577. Horizontal bars mark the outer horizon, where V — 1. Light from the outside universe is 
visible only from outside the Cauchy horizon, so lines terminate infmitesimally outside the Cauchy horizon even in models in 
which the baryons drop inside the Cauchy horizon (middle model). 



black hole, shown in the left panel of Figure ^1 is 
similar to that of the corresponding vacuum black hole, 
the Schwai'zschild solution. As the observer plunges to 
the central singularity at zero radius, the black hole 
increases monotonically in angular size, reaching 90° at 
the singularity. The blueshift of photons at the edge 
of the black hole tends to infinity. The appearance 
can be attributed to the enormous tidal force near 
the singularity. The same tidal force that stretches 
the infaller radially and crushes them horizontally also 
aberrates photons so that they appear to come from 
a thin horizontal blueshifted band of light on the 
sky. Photons from above the observer are redshiftcd, 
while photons from the thin horizontal band are highly 
blueshifted. Although the blueshift at the edge of the 
black hole tends to infinity as the observer approaches 
the singularity, the amount of time that the observer sees 
pass by in the outside universe, the integral of blueshift 
over proper time, is finite. 

The appearance of a charged, non-conducting black 
hole, shown in the middle panel of Figure 1121 is again 
similar to that of its vacuum counterpart, the Reissner- 
Nordstrom (RN) solution, down to the inner horizon. 
From the point of view of an outgoing observer, such as 
one in the rest frame of the charged baryonic fluid, the 
black hole increases in angular size until, at the Cauchy 
horizon, the black hole covers the entire sky. The view of 



the outside universe correspondingly decreases in size and 
becomes brighter and more blueshifted, until the view 
disappears at the Cauchy horizon in an infinitely bright, 
blueshifted, concentrated flash. The entire future of 
the outside universe passes in that infinitely blueshifted 
flash. It should be commented that such an illuminating 
experience could not actually happen in reality, because 
photons from the outside universe are necessarily ingoing, 
and the presence of even the tiniest amount of ingoing 
matter or radiation will prevent passage through the 
Cauchy horizon — see §V A of Paper 2. 

From the point of view of an ingoing observer on the 
other hand, such as one who free-falls radially from zero 
velocity at infinity, the charged, non-conducting black 
hole (middle panel of Figure ^1 again) first increases 
in angular size, then shrinks as the observer approaches 
the inner horizon. Just above the inner horizon, the 
angular size of the black hole is finite, and the blueshift 
is still finite, just as it is in the RN geometry. In the 
RN solution, an ingoing observer would pass through the 
inner horizon, and in the instant of passage would see 
a new image of the outside universe, reflected by the 
gravitationally repulsive core of the black hole, appear 
at the center of the black hole in an infinitely blueshifted 
flash containing the entire history of the universe. In the 
similarity solution by contrast, the ingoing observer does 
not pass through the inner horizon, and no new image 
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of the universe appears in an infinitely blueshifted flash. 
Instead, the ingoing observer falls to zero radius while 
still remaining just outside the inner horizon. As the 
ingoing observer falls to smaller radius, they encounter 
outgoing baryons accreted at ever earlier times, back 
to the initial collapse event at r = 0. In this regime 
where the ingoing observer remains close to the inner 
horizon, the appearance of the outside universe evolves 
simply by becoming more and more Lorentz-boosted. 
Although the blueshift tends to infinity as the ingoing 
observer approaches zero radius, the amount of time that 
the ingoing observer sees pass by in the outside universe 
remains finite. This is unlike the outgoing observer, who 
(if they could actually see it) would see the entire future 
of the universe pass by as they cross the Cauchy horizon. 

The right panel of Figure El shows the appearance 
of the discretely self-similar charged baryonic model, 
with conductivity equal to the continuum threshold 
conductivity k = Koo , from Figure If 01 As is evident 
from Figure El the oscillations in the model are mildy 
visible to a free-fall observer, but, curiously, practically 
invisible to an observer at rest in the baryonic frame. 
In the baryonic frame, the apparent angular radius of 
the black hole remains near 90° all the way from the 
outer sonic point to the singularity at zero radius. Near 
the singularity, the 90° size can be attributed to the 
enormous tidal force, as in the Schwarzschild geometry, 
but at large radii the near 90° size is presumably 
fortuitous, a happenstance of the relativistic aberration 
from the motion of baryons at and inward of the sonic 
point. As in the case of the uncharged black hole, 
although the blueshift at the edge of the black hole tends 
to infinity as the observer approaches the singularity, the 
amount of time that the observer sees pass by in the 
outside universe remains finite. 

For those interested in carrying out approximate ray- 
tracing inside one of the black holes described in this 
paper, it is useful to note that the function H(X), which 
plays the essential part in ray-tracing, equation f7U|) . is 
reasonably approximated as a cubic polynomial in X if 
the baryons collapse straight to a central singularity, or 
as a quartic polynomial in X if the baryons drop through 
the Cauchy horizon. The approximation becomes exact 
in the asymptotic limit of small accretion rates and small 
conductivities, where H(X), equation |75]l. becomes 
a cubic polynomial if the black hole charge is zero 
(Schwarzschild limit), or a quartic polynomial if the 
black hole charge is non-zero (Reissner-Nordstrom limit). 
Usefully, the approximation of H(X) as a cubic or quartic 
polynomial is not bad even at finite accretion rates and 
conductivities. 



VI. SUMMARY 

In this paper we have investigated self-similar solutions 
for spherically symmetric charged black holes that 
accrete a relativistic fluid (pb/pb = 1/3) of charged, 



electrically conducting baryons. The solutions are 
characterized by three free parameters: the accretion 
rate, the charge-to-mass of the black hole (generated 
self-consistently by the charge-to-mass ratio of accreted 
baryons), and the electrical conductivity. We do 
not require that solutions be regular at zero radius, 
but rather integrate inwards from boundary conditions 
established at a regular sonic point outside the outer 
horizon of the black hole. 

The accreted charged baryons undergo one of two 
possible fates: either they plunge directly to a spacelike 
singularity at zero radius, or else they drop through a 
Cauchy horizon. In the latter case the baryons probably 
undergo a shock just inside the Cauchy horizon, but, 
whether or not a shock occurs, the similarity solutions 
do not continue consistently to zero radius, but rather 
terminate at an irregular sonic point where the proper 
acceleration diverges. We argued in WV Dl that the failure 
of the similarity solutions to continue to zero radius inside 
the Cauchy horizon does not invalidate the solutions as 
a whole, because the failure is hidden behind the Cauchy 
horizon, and cannot be communicated to the solution 
outside the Cauchy horizon. It remains a matter for 
future numerical investigation to discover what really 
happens beyond the Cauchy horizon. 

In the solutions where the baryons drop through the 
Cauchy horizon, the geometry inside the Cauchy horizon 
differs in one crucial respect from the corresponding vac- 
uum solution, the Reissner-Nordstrom (RN) geometry. 
Whereas in the RN geometry the radial 4-gradient is 
spacelike inside the Cauchy horizon, in the similarity 
solutions the radial 4-gradient is timelike. This means 
that in the similarity solution it is impossible to remain 
at rest (at fixed radius r) inside the Cauchy horizon: 
all locally inertial frames necessarily fall to smaller 
radii. Because the similarity solutions do not continue 
consistently to zero radius, we cannot say exactly what 
happens, but if the radial gradient remains spacelike all 
the way to zero radius, then the baryons will fall to zero 
radius, where there is presumably a spacelike singularity. 

It may seem paradoxical that the spacetime inside 
the Cauchy horizon is subluminal, yet nevertheless the 
baryonic fluid must flow to smaller radii. The resolution 
of this paradox is that from the point of view of 
observers inside the Cauchy horizon, the Cauchy horizon 
is actually contracting to smaller radii. Notwithstanding 
this contraction, observers who fall through the Cauchy 
horizon at later times find the Cauchy horizon at larger 
radii. A consistent way to understand this is to think 
of coordinate time as running backwards inside the 
Cauchy horizon (though of course proper time always 
runs forward in the usual way). The topsy-turvy world 
inside the Cauchy horizon was discussed in fliVBI 

The question of whether the baryons plunge to a 
singularity or drop through the Cauchy horizon depends 
on the parameters, as summarized immediately below. 

If the black hole is uncharged (because the baryons 
it accretes are uncharged), then the baryons plunge to 



20 



a central spacelike singularity, as in the Schwarzschild 
solution. 

If the black hole is charged and the conductivity 
is zero, then the accreted baryons, repelled by the 
charge of the black hole generated by previously accreted 
baryons, naturally become outgoing, and drop through 
the outgoing inner horizon, the Cauchy horizon. 

If the black hole is charged and if the dimensionless 
electrically conductivity k exceeds a certain threshold 
value (whose value depends on the accretion rate 
and charge-to- mass ratio), then the baryonic fluid 
can neutralize inside the black hole, and plunges to 
a spacelike singularity. If on the other hand the 
conductivity is less than the threshold value Koo, then 
in most cases the baryons drop through the Cauchy 
horizon. However, there is also a discrete spectrum 
Kn (n = 1,2,...) of conductivities, asymptoting to the 
threshold k, OCi at which the baryons plunge to a central 
singularity rather than dropping through the Cauchy 
horizon. Solutions between the discrete values Krfi are 
characterized by the number of times that the baryonic 
fluid oscillates between ingoing and outgoing before 
dropping through the Cauchy horizon. 

At the threshold conductivity n^, the baryonic fluid 
oscillates between ingoing and outgoing an infinite 
number of times before hitting the central singularity, 
and the solution exhibits a discrete self-similar behavior 



reminiscent of that observed in critical collapse [2^, |2^, 

m. 

In <JV] we discussed what an observer who falls inside 
one of the black holes considered in this paper would see. 

Mass inflation at the Cauchy horizon [27I l66| does 
not occur in the solutions considered in this paper. 
Mass inflation requires the presence of both ingoing 
and outgoing fluids near the inner horizon, and in this 
paper we have deliberately chosen to consider only a 
single baryonic fluid. Although a finite conductivity does 
allow oppositely charged fluids to diffuse through each 
other, this is not enough to permit mass inflation. In a 
companion paper, Paper 2 [33, we allow the black hole 
to accrete not only baryons but also dark matter, and 
it will be seen that streaming between baryons and dark 
matter leads to mass inflation. 
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APPENDIX A: GENERAL SPHERICALLY SYMMETRIC TETRAD 



This appendix presents expressions for the vierbein, metric, connection coefficients, and Einstein tensor, for the 
most general spherically symmetric orthonormal tetrad. 

It is convenient to proceed in two steps, carried out in S ^A II and IA 21 immediately below. The first step is to 
generalize the vierbein adopted in the text, equations I©, to the case of arbitrary radial coordinate. The second step 
is to complete the generalization. 



1. Arbitrary radial coordinate 

In the text, the gauge choices v = fj, = and A = 1, equations (JSJ), were imposed on the general vierbein 10}, 
bringing the vierbein to the form JfjJ) . The gauge choice A = 1 is equivalent to setting the radial coordinate r to be the 
circumferential radius, defined such that the proper circumference of a circle is 27rr. In this subsection of Appendix lAl 
the gauge choices v = fx = are retained, but A is permitted to be arbitrary, so that the radial coordinate is arbitrary. 

Denote the arbitrary radial coordinate r 1 and the corresponding vierbein coefficients a' , /?', 7', and A' with primes, 
to distinguish them from the circumferential radial coordinate r and vierbein coefficients a, (3, 7 considered in the 
main text. The vierbein e™^ are then 



eo = a , a = , e z = $ Xi , e t J = 7' x&j + A' - XiXj) (Al) 



with inverse 



e% — — 7 , e°i — , e i Q — - — - Xi , e? i — — XiXj + — (Sij — XiXj) . (A2) 
a' a'y 7' A' 



The metric is 



, 2 dt 2 1 / , p'dt\ 2 r' 2 do 2 

ds 2 = - + —[dr'- — — + 5- . A3) 

a' 2 i 2 \ a J A 



It is evident from the metric (IA3|) that the quantity r defined by 
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is the circumferential radius r of the main text. A relation between the new vierbein coefficients a', (3', 7', A' of 
equation (|A1() and the vierbein coefficients a, f3, 7 in the main text is most easily established by noting that the 
directed derivatives d mi equation (JHJ, are independent of the choice of coordinate system x M , and in particular are 
independent of the choice of radial coordinate. Thus the directed derivatives in the time t and radial r directions 
must be the same irrespective of the choice of radial coordinate, r or r'\ 



1 9 
dt=a dl 

fl 9 



■a 



d_ 

dr 



«' ' 9 
9 t = a dt 



Q/ / 9 



(A5) 



(A6) 



It follows from equalities (|A5|) that a, /?, 7 in the main text are related to a', f3' , 7', A' by 

The inverse formulae for j3' and 7' in terms of (3 and 7 are /?' = /3A' + r 5 t A' and 7' = 7 A' + r d r \' . All the formulae 
given in the main text now carry through with r, a, f3, and 7 given by equations l|A4|l and l|A6[) . One useful relation to 
note is that the proper acceleration g' and proper velocity gradient h! defined analogously to g and h of equations 121|) 
and J22J), 



d r In a' , h! 



9 P' o , , 

57 ln7 



/3' 



, 9 In a' 
<9r' 



(A7) 



are the same as their unprimed counterparts 



.9 = .9 



(A8) 



2. General case 

The most general form of the vierbein e m ' 1 consistent with spherical symmetry is given |69j by equation The 
general vierbein can be obtained from the vierbein of the previous subsection, equation IjAljl . by boosting the tetrad 
frame radially by boost angle £, and rotating the tetrad frame about the radial direction by angle £. The boosted and 
rotated vierbein coefficients (doubly primed, to distinguish them from the singly primed coefficients of the previous 
subsection) are then 

a" = a' cosh f , v" = a' sinh f , /?" = /3' cosh £ + 7' sinh f , 7" = 7' cosh £ + /3' sinh £ , A" = A' cos (, fj," = A' sin ( . 

(A9) 

Equations l|A9() imply the relations 



/2 //2 //2 /2 „/2 i/2 „//2 / / // n all II I pi II all II II \/ 2 \// 2 i " 2 / A 1 n\ 

a = a — v , 7 — p =7 — p , 0:7 = 07— p^ , ap=ap— 7^ , A = A +/i (All)) 

which are sufficient to determine the singly primed, unboosted, unrotated vierbein coefficients in terms of the doubly 
primed, boosted, rotated vierbein coefficients. The vierbein e m ^ are [the following is the same as equation Jl}, but 
with doubly primed coefficients] 

e ° = a" , e,° = v" x< , e * = f3" £< , 

ei 3 = 7" sc^ + A" (<5.y - XiXj) + n" SijkXk (AH) 



with inverse e m M 



V ;/' Pi" 

_ 7 ~ f i _ P 



a"7" - /3'V ' a"7" - /3'V ' a"j" - f3"v" 

3 - °" - - A " /i" 

6 1 _ a"7" - [3"v" XiXj + A" 2 + /u" 2 ~~ ^ ~ A" 2 + /x" 2 ' 

Boosting and rotating the tetrad has no effect on the coordinates x^, and therefore no effect on the metric, which 
continues to be given in terms of a', f3' , 7', and A' by equation i|A3l) . In terms of the boosted and rotated vierbein 
coefficients a" , v" : (3" , 7", A", and //', the metric is 



<ll- a" 2 -v" 2 r., K/3" - yv)<ft 



a" 2 - //" 2 (a'V - /?'V) 2 



dr' 



a" 2 - i/" 2 



'do 3 



" 2 -L- //" 2 



A" z + M' 



(A13) 
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with determinant 



(a'V - P"u") (A" 2 + /j," 2 ) 



(A14) 



The directed derivatives d' t ' and d" in the time i and radial r directions are Lorentz boosted versions of the directed 
derivatives dt and d r of the previous section, equations (|A5(1 : 



d' t ' 
d'l 



a"^-+/3"^7 = cosh£d t +sinh£d r 



at 
o 



dr 

a 



7"ttt — smn £ dt + cosh £ 9,- 



dt 1 dr 

In terms of doubly primed vierbein coefficients, the tetrad connection coefficients are explicitly 
~dv" di' 



L/u L/^> „ i u II rill /Is 

~~dT ~dr r ~ 4 ^ ~ ^ ) 



r 



^-iarin(A" 2 +/i" 2 ) 



f A M a// i 



A" z 



A" 



Is; In (A" 



2 ' M " 2 



„ . AV 
XkXl A" 2 + i 



9" In 



,2 r 



+ (4i - X k Xl) —7 



(A15a) 
(A15b) 

(A16a) 

(A16b) 

(A16c) 
A*' 



(A16d) 



Tensor quanties follow straightforwardly by boosting and rotating the tetrad frame from the v = 0, fi = frame of 
the previous subsection, ijA II For example, the Einstein tensor is 



G mn = 2[(R+ S) t m t n + F (i m f n + f m t n ) + P f m f n + S (if 



(A17) 



where F, R, P, S are given by equations (|25|l [with r, a, (3, and 7 being unprimed quantities, given in terms of 
primed quantities by equations (|A4|) and l|A6(l . and thence in terms of doubly primed quantities by the inverse of 
equations (|A9|l ] , and the 4- vectors t m and f™ are respectively the unit time vector and the unit radial vector boosted 
by boost angle £ out of the v = 0, /1 = frame 



r = cosh£ 



sinh£ 



sinh £ cc, , 
= cosh £ Xi . 



(A18) 
(A19) 



Note that the rotation of the tetrad frame by angle C about the radial direction leaves the Einstein tensor G mn , and 
the energy- momentum tensor T mn , and more generally any symmetric tensor of rank 2, unchanged. 



APPENDIX B: PROOF THAT j 2 - (3 2 > AT THE OUTER SONIC POINT 



This Appendix proves, rather unprettily, the assertion 
in fllll HI that, for the sonic point at the outer boundary 
to be causally connected to (not separated by a horizon 
from) an asymptotically flat empty region of space at 
large radius r, it is necessary that the radial 4-gradient 
(/3, 7) be spacelike at the sonic point, 7 2 — 1 > 0. The 
proof applies to the self-similar solutions not only of the 
present paper, but also of Paper 2, where the black hole 
accretes dark matter in addition to baryons. Quantities 
below subscripted d refer to dark matter. 

The assertion will be proven for the case where the 
geometry consists of a spherically symetric self-similar 
region surrounded by empty space. As remarked in 



mil HI the self-similar solutions modeled in this paper 
do not necessarily continue self-consistently to infinite 
radius. We suppose therefore that the self-similar region 
is truncated at some point outside the outer sonic point, 
and that space beyond this truncation point is empty. 
Actually the proof works if the region outside the self- 
similar region is not empty; it suffices to assume that 
7 2 — 1 > at the outermost point to which the the 
self-similar solution is considered to extend. 

According to the definition of interior mass, 7 2 — 1 — 
1 — 2M/r, equation (|24p. In the hypothetical empty 
portion of space at large r, the geometry is necessarily 
Reissner-Nordstrom. By assumption, the empty region 
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is outside any horizon, which implies that j 2 — f3 2 > in 
the empty region. 

Now consider the self-similar region of space. The 
differential equation governing the self-similar evolution 
of 7 2 — 1 is, in the baryonic frame, 

17 J ' = 2T [(1+^)^7 + ztftld] (Bl) 

which includes a contribution from pressureless dark 
matter (Paper 2). To prove that 7 2 — (3 2 remains 
positive, assume the opposite, that there is a point 
where r ) 2 — 1 drops through zero from positive to 
negative. At the point where j 2 — /3 2 passes through zero, 
the derivative on the right hand side of equation (|B1|) 
must be negative, so consider the various terms in this 
expression. As remarked in flIII CI consistency of the 
self-similar solutions requires that £ r always remain of 
the same positive sign, since if £ r changed sign it would 
indicate that the baryonic fluid turns back on itself, 
which physically cannot happen. The dimensionless 
proper baryonic and dark matter densities z and z<j in 



equation (|B1|) are necessarily positive. By assumption, 
the dark matter is causally connected to, and indeed falls 
from, the outer empty region of space, where jd > 0. 
Since ^ 2 — (3 2 — ^ 2 — (3 2 , it follows that j d > throughout 
the region where j 2 — /3 2 > 0. Since the baryonic 4- vector 
((9,7) is related to the dark matter 4-vector (/3d, 7d) by 
a Lorentz transformation, equation (5) of Paper 2, it 
follows likewise that 7 > throughout the region where 
7 2 — (3 2 > 0. For the derivative on the right hand side of 
equation (|B1|) to be negative thus requires that at least 
one of (3 or (3d be negative. Suppose without loss of 
generality that it is (3 that is negative. If /3 < 0, then 
7-/3 > 7 + (3, and for 7 2 — /3 2 to pass through zero, it 
is necessary that 7 + j3 = 0. Now 7 + f3 = 7(1 — V) + rj 
according to equation for the velocity V. But 7 > 0, 
and V < 1 since the point is outside the horizon, so 
7 + (3 > i]. But i] — l/£ is a strictly positive quantity 
everywhere outside the horizon. Thus 7 + [3 is strictly 
positive, contradicting the proposition that 7 + (3 = 0. 
This proves the theorem. 



